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Theorie dreidimensionaler Setzungsvorginge in Tonschichten 


Von G. Heinrich und K. Desoyer 


1. Einleitung. In der bisherigen Literatur wurden die Konsolidierungs- und Setzungsvorgange 
_in Tonschichten in der Regel nur eindimensional behandelt!2. Die wenigen Ansatze zu einer 
_ mehrdimensionalen Setzungstheorie versuchen meist, die Porenwasserstrémung und die Setzung 
_ der Tonschichten voneinander unabhingig zu erfassen®. Dies liefert zufolge der Koppelung dieser 

beiden Vorginge keine verlafliche Grundlage. In einer Reihe von Arheiten‘ gibt Biot eine Er- 
_ weilerung der Setzungstheorie auf mehrdimensionale Vorgange. Da er jedoch nicht von einer 
_ Krafteanalyse fiir die beiden Komponenten (Festkérper, Fliissigkeit) und den beiden zugehérigen 
_ Kontinuitatsgleichungen ausgeht, erscheinen seine Ansitze nicht aus den Grundlagen entwickelt. 
_ Er gelangt zu dhnlichen Ausgangsgleichungen wie in der vorliegenden Arbeit, behandelt aber, 
_ soweit den Verfassern bekannt ist, damit kein konkretes raéumliches Setzungsproblem. 


_ 2. Problemstellung, Voraussetzungen und Ausgangsgleichungen. Der Gleichgewichtszustand 
eines wasserhaltigen Bodens wird gestért. Der durch diese Stérung hervorgerufene Verlauf der 
_ Grundwasserstrémung und der Bodenbewegung soll ermittelt werden. Dabei wird vorausgesetzt: 


1. die Poren des Bodens sind stets vollstandig mit Wasser gefiillt; 
2. Wasser und Festkérpermaterial sind inkompressibel; 


: 3. die zwischen Wasser und Boden auftretenden Molekularkrafte werden nur durch Korrektur 
_ des Porenvolumens beriicksichtigt’; 


4. die Filterwirkung des Festkorpers sei geniigend grofs, so daB eine ,,quasistationare Bewegung“, 
bei der die Tragheitskrafte in der Flissigkeit und im Festkérper gegeniiber den anderen Kraften 
_ vernachlassigt werden kénnen, bereits nach einem verschwindend kleinen Zeitintervall nach Be- 
- ginn der Stérung auftritt; 


5. es bestehe ein linearer Zusammenhang zwischen dem Vektor der Filtergeschwindigkeit und 
_ dem Druckgradienten in der Fliissigkeit im Sinne des Darcy-Gesetzes; 


6. fiir nicht zu groBe zusatzliche Deformationen des Festkérpers bestehe nur eine lineare Ab- 

_ hangigkeit zwischen der Anderung des Deformationstensors und der Anderung des Spannungs- 

_ tensors des Festkérperskelettes, solange die zeitliche Dilatationsinderung ihr Vorzeichen nicht 
wechselt. . 


Die im folgenden verwendeten Begriffe Festkérpergeschwindigkeit (Skelettgeschwindigkeit), 
_ Filtergeschwindigkeit, relatives Porenvolumen, Korn-zu-Korn-Krafte, Druck in der Flissigkeit 
Bice: sind im Sinne eines Makromodells aufzufassen, das durch Mittelwertbildungen aus den tat- 
- sachlich auftretenden, aber im einzelnen nicht bekannten Mikrowerten hervorgeht® ’. 

- Yur Ableitung der Kontinuitatsgleichungen fiihren wir folgende Groen ein: n das auf die 
 Volumeinheit des Gemisches (Wasser und Festkérper) bezogene Porenvolumen, Ux der Geschwin- 
igkeitsvektor des Festkérpers, » der Vektor der Filtergeschwindigkeit relativ zum Festkérper, 

lp der Vektor eines Oberflaichenelementes. 


1 Terzaghi-Jelinek, Theoretische Bodenmechanik, Berlin 1954. : 
_ 2 Terzaghi-Fréhlich, Theorie der Setzung von Tonschichten, Leipzig und Wien 1936. 
3 N, Carrillo, J. Math. Phys. 21 (1942) S. 1. : ; 

4 M.A. Biot, J. Appl. Phys. 12 (1941) S. 155 und S. 426; M. A. Biot u. F. H. Clingan, J. Appl. Phys. 12 
(1941) S.578; M. A. Biot, J. Appl. Phys. 26 (1955) S. 182; J. Appl. Mech. 23 1, (1956) S. 91; M. A. Biot u. D. 
_ G. Willis, J. Appl. Mech. 24 (1957) S. 594. : : 
5 J, Kozeny, Hydraulik, S. 381, Wien 1953. __ 

G. Heinrich u. K. Desoyer, Ing.-Arch. 23 (1955) S. 73. 
7 G. Heinrich u. K. Desoyer, Ing.-Arch. 24 (1956) S. 81. 
Ss ee ee cone 16 
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Wir denken uns eine raumfeste Hiille 0, die das Volumen V umschlieBt, gemaB Abb. 1. Das 
Volumen des im Zzitelement dt durch diese raumfeste Hille hindurchtretenden Festkérpermaterials 
ist!, [(l—n)b,-dodt, die Abnahme des Festkérpervolumens innerhalb der Hille ist 

é 


—= [a —n)dV dt. Wegen Voraussetzung 2. gilt daher 
V 


[Q—n) 0g do de =—F, (Un) ar at 
) 4 


oder nach Anwendung des Gaufschen Satzes 


[71a —») ol — Gav =o. 
Abb. 1. Zur Ableitung V 
der Kontinuitatsgleichungen. c 
Da diese Gleichung fiir jedes beliebige Volumen 
gelten mu, folgt daraus 
én 


(l—n)V-bg=24(Vn)-de. — () 


Abb. 2 (rechts). Zur Definition 
der Standrohrspiegelhéhe h. 


Um die Kontinuitatsgleichung fiir die Flissigkeit zu erhalten, denken wir uns wieder eine 
raumfeste Hiille (Abb. 1). Wiirde das Wasser relativ zum Festkérper in Ruhe bleiben, so ware 
nx: do dt das in der Zeit dt durch das raumfeste Oberflachenelement do transportierte Flissig- 
keitsvolumen. Besitzt die Fliissigkeit noch zusatzlich eine Relativbewegung gegen den Festkérper, 
so vermehrt sich das in der Zeit dt durch do hindurchtretende Flissigkeitsvolumen um b -: dp dt. 

Es gilt daher gema8 Voraussetzung 2. die Beziehung 


7) 
[ere +») -ao dt =—,, \ndVd. 
6 Vv 
Daraus folgt nach Anwendung des Gaufschen Satzes 


V-(nbdx +b) +2 =0 


oder 
nV te + 0-0 =—[F + (n)- Og) (2) 
Aus (1) und (2) ergibt sich 
V-(be +b) =0. (3) 
Fiir die relative Filtergeschwindigkeit » gelte das Darcy-Gesetz in der Form 
b=—kTh. (4) 
Hierin ist k die Durchlassigkeit und h die Standrohrspiegelhéhe gema® der Beziehung 
U 
bee (5) 


U ist hierin das Potential der Massenkraft (Schwerkraft), bezogen auf die Masseneinheit, p der 
statische Druck in der Flissigkeit, yy das spezifische Gewicht der Fliissigkeit, g die Fallbeschleu- 
nigung. 

Mit den Zahlrichtungen nach Abb. 2 ist U = — gz und daher 


1 
Mts Ys BA) = Pl, ¥, 21) — 2 - (5a) 


Wir denken uns nun durch cine kiérperfeste Hiillflache ein endliches Volumen abecoreac 
Dann gilt die folgende Gleichgewichtsbedingung fiir den Festkérperanteil2 


VG =AVp +8 —*Ue— Ww) 7y + pyA. (6) 


1 Der Punkt in Zeilenmitte bedeutet das skalare Produkt i P i i i i 
peteiahdar dneeerile Peas , zwei Punkte in Zeilenmitte nebeneinander 


2 Siehe Gl. (9) der in FuBnote 7 von Seite 225 zitierten Arbeit. 
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Die Hiullflache sei dabei so gelegt, daB kein Korn des Festkérpermaterials durchschnitten 
wird. In diesem Fall bedeutet © den Spannungsaffinor der Korn-zu-Korn-Krafte; 7 ist fiir ein 
makroskopisches Flaichenelement der Hiillflache das Verhiltnis des von Flissigkeit bedeckten 
Anteils der Flache zur gesamten Flache (A = 1 fiir Punktberiihrung der Kérner). yx ist das spe- 
zifische Gewicht des Kornmaterials. Die obige Schnittfiihrung ist dadurch begriindet, daB laut 
Voraussetzung 2. die Deformation des Festkérpers nur durch Anderung des Porenvolumens ver- 
ursacht wird und diese nur auf das Wirken der Korn-zu-Korn-Kriifte zurtickgefiihrt werden kann. 
Wiurde man auch elastische oder plastische Deformationen der Kérner in Betracht ziehen, dann 
. hierfiir' der durch einen statistischen Schnitt gewonnene Spannungstensor verwendet 
werden. 

Laut Voraussetzung 6 gelte folgender Zusammenhang fiir die durch die Stérung verursachten 
zusatzlichen substantiellen Anderungen 2’ des Deformationstensors® des Festkérperskelettes 
und © des Spannungstensors © der Korn-zu-Korn-Kriafte: 

Y =a[(l +7) 6’ —» SG]. (7) 
Darin bedeuten « und v zwei das Deformationsyerhalten des Skelettes charakterisierende positive 
Konstante? (in Analogie zum reziproken Elastizitatsmodul und zur Querzahl), © den Einheits- 
affinor und S{ 
Ss 6’ ie & b) 
die erste Invariante von G’. 

Man kénnte Gleichung (7) noch durch Einbeziehung des Affinors der Deformationsgeschwin- 
digkeiten erweitern. Wir wollen hier aber bei den bisher tiblicherweise angenommenen physi- 
kalischen Zusammenhangen bleiben und diese nur in geometrischer Hinsicht auf dreidimensionale 
Falle verallgemeinern. 

Aus (7) 1aé8t sich auch die erste Invariante 


Di =e’ (Volumsdilatation) (8) 
des Verzerrungsaffinors D’ mit S; in Zusammenhang bringen. Man erhilt 
é-= DE =D -- © =a —27) S,. (9) 


In sinngemaBer Erweiterung der bisherigen Ansatze? wird angenommen, daf die positiven 
Beiwerte « und y nur von sgn De’/dt abhingen, wenn D/dt die substantielle Ableitung beziiglich 


der Skelettbewegung bedeutet. Day < = sein muB, gilt sgn De’ /dt = sgn DS; /dt. Ist De’ /dt negativ 


(Konsolidierung), dann sollen die Beiwerte mit «, und »,; ist De’/dt positiv (Schwellung) mit a, 
und v, bezeichnet werden. Es ist anzunehmen, da die Beiwerte « und v in geringem MaBe auch 
noch in anderer Weise von den Invarianten des Deformationsaffinors und des Affinors der De- 
formationsgeschwindigkeiten abhangen werden, doch sollen im folgenden «,, »,, «,, 7, als Kon- 
stante angesehen werden. Ferner soll angenommen werden, da der Sprung der Wertea und ¥ 
nur bei Vorzeichenanderung von De’/dt auftritt. 

_Aus (7) folgt mit Verwendung von (9) 


6 : 


v oes 
Fe pera (10) 
Um einen Zusammenhang zwischen der Anderung n’ des relativen Porenvolumens n und der 
Dilatationsainderung ¢’ herzustellen, denken wir uns ein Volumelement dV, das sich aus den An- 

teilen dVx (Kornmaterial) und dV, (Wasser) zusammensetzt: 


Nach Definition ist 
rs Vy 12 
= rie (12) 
Bedeutet D die substantielle Anderung, dann ist 
_ DV) a 
De rt dV. ° ( ) 
~ 1 Siehe FuBnote 7 von Seite 225. 
2 Siehe Gleichung (63). 
3 Siehe FuGBnote 1 und 2 von Seite 225. 
16* 
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Wegen der Erhaltung des: Kornvolumens folgt aus (11) 
D(dV) = D(dV yp) (14) 
und damit aus (12) 
Dn: dV +n D(dV) = DV yy) . 
Daraus ergibt sich bei Verwendung von (13) und (14) 
Dn 


1l—n 


en 


Die Integration liefert, wenn zum Anfangswert éy das relative Porenvolumen ny gehort, 


ie Comers —=E—£- 


1— nm 


Setzt man 
ae) (18) 


so ergibt sich 


In (1 =z = et (16) 


Da die Gleichungen (7) und (9) nur fiir kleine Anderungen Giiltigkeit haben, wird man im 
Rahmen dieser GréBenordnung statt (16) die Naherung 
n’ , 
ees (17) 


verwenden. 


3. Linearisierung der Ausgangsgleichungen. Die Gleichungen (1), (3), (4), (5), (6), (7) und (17) 
sind die Ausgangsgleichungen des Problems. | 

Wir gehen nun von einem beliebigen Anfangszustand aus. Die zugehérigen GroBen seien mit 
dem Index 0 bezeichnet. Zur Zeit t = 0 werde eine ,,kleine Stérung“‘ eingeleitet, die in dem Sinn 
als klein anzusehen ist, daB die Produkte der StérgréBen gegenitiber diesen selbst vernachlassigt 
werden kénnen. Die Stérgré8en sollen durch einen Strich gekennzeichnet werden. Auf Grund 
der obigen Annahme kénnen nun die Ausgangsgleichungen linearisiert werden. Die Gleichungen 


(1), (3), (4), (5) und (6) gehen dann iiber in 


—n' V+ bg.) + (1m) (V+ 0g) = + (Vg: Vg) + Vn’ + dx.) 5 (la) 
V-(0, +b’) =0, (3a) 
py’ = — (ky Vio EVI) (4a) 
he 
War, (5a) 


VG = AyVp' +i Vey—"S8—™O VU + p' Why + PoVH. (6a) 


Dazu kommt noch die unveranderte Gleichung (7). 

Ist der Ausgangszustand ein Gleichgewichtszustand ohne Grundwasserstrémung, sind ko, ny 
und A, nicht vom Ort abhangig und ist iiberdies 2’ vernachlassigbar klein! gegen J), dann verein- 
fachen sich die Gleichungen (la), (4a) und (6a) und man erhalt das folgende Gleichungssystem: 


(l—m) (7-0) ==, (1b) 
V «(be +’) =0, (3b) 
b’ = —k,Vh’, (4b) 
gene 
yw (ch) 


1 Diese Annahme ist dadurch gerechtfertigt, da bei Tonen nahezu Punktberiihrune der Ké - 
handen ist, so daB A nicht wesentlich vom Wert 1 abweichen kann. " nal sae ie ma Aue 
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VS =), 7p’ —" Ve tw) py, (6b) 
Y all +) 6 SG], (7b) 
— ae (17b) 


Dies sind 16 skalare Gleichungen, die zusammen mit vorgegebenen Anfangs- und Randbedin- 
gungen das Problem vollstandig bestimmen. Man kann nadmlich die sechs Komponenten von 9’ 
durch die drei Verschiebungskomponenten &’,n’,¢’ des Festkérpers in bezug auf ein raumfestes 
Koordinatensystem (x, y,z) ausdriicken. (Lagrangesche Darstellung.) Die Bewegung der Fliissig- 
keit relativ zum Festkérper kann dann in Eulerscher Darstellung durch die Komponenten v,, Us» 
v, der Filtergeschwindigkeit beschrieben werden, Als die 16 unbekannten Funktionen kénnen 
also die 10 GréBen &, 7’, 2’, v., U,, V,, 0’, &, p’, h’ und die sechs Komponenten von 6’ angesehen 
werden. 


4, Umformungen des gewonnenen Gleichungssystems. Es soll zundchst eine Differentialglei- 
chung fiir e’ abgeleitet werden. Aus (7b) folgt 


V-D’ =o(1+7)V: SG —arVS, 
oder mit Verwendung von (6b), (9) und (17b) 


V-D =a(1 +7) Yo 1 rete SAA ( eae) eae (18) 


Daraus erhalt man mit 4U = 0 


V-V-9) =a0(1 +») 0 Ap’ — "S28 (1 — 1m) (VU): (Ve) — 


v , 

Toole * (19) 
Bedeutet u’ den Verschiebungsvektor des Festkérpers an einer Stelle mit dem Ortsvektor r, 

dann gilt! 


Y (ce tu)y=>5C0.w +u47). (20) 


Weiter ist 
D (vctu)=D(r) +u-V,D) +:°°. 

Wir setzen nun wie ublich im Sinne einer Theorie erster Ordnung voraus, daf die Verschie- 
bungen klein sind gegen die Abmessungen des Bereichs. Dann kann man naherungsweise D (x ere u’) 
durch D’(r) ersetzen und alle Beziehungen durch die urspriinglichen unabhangigen Variablen 
ausdriicken. Die substantielle Ableitung in bezug auf die Skeletthewegung kann bei dieser Naherung 
durch die partielle Ableitung nach der Zeit ersetzt werden. Dann geht (20) tiber in 


DS = Viw +uV), (21) 
und daraus folgt 
VO = 5 [aw +00 -v). (21a) 
Mit 
; Wise (22) 
kommt 


VV “d= Ae. 
Damit und mit (5b) erhalt man aus (19) 


Ae = OFC) Th yy AW’ — "R= EE (1 — my) (VU) + Ve’) « (23) 
Aus (4b), (3b), (1b) und (17b) folgt 
Ae (24) 
ky Ot 
Aus (23) ergibt sich mit (24) schlieBlich 
ee (25) 
Ae’ +%,(VU): Ve’) =x ar? 


1 Komma bedeutet das dyadische Produkt. 
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mit den Abkirzungen 


__ a(l1+y)(1— 2») (1 —™) ¥e—vm) (26) 
ire <9 Seay g ; 
4 (1+ 7%) — 2) dove 27 
See ae ee (27) 
Wegen (9) gilt auch we 
AS, +4 (VU): (VSi) =x. see 


Um eine Differentialgleichung fiir h’ allein zu bekommen, kann man Gleichung (25) partiell 
nach der Zeit differentiieren und erhalt dann mit (24) 


A Ak’ +2 (VU) + [(V(Ah’)] =x = (AN). (29) 
Die GréBe Ah’ gehorcht demnach derselben Differentialgleichung wie ¢’. 


Zu einer Beziehung zwischen uw’, e’ und h’ gelangt man durch Gleichsetzen der rechten Seiten 
von (18) und (21). Mit (22) und (5b) erhalt man bei Verwendung der Abkiirzungen (26) und (27) 


/ 2(— Y) 1 / / e-2 , 
AN = ass; 2a—n le | 4,6’ VU —kyxVh'| . (30) 


Man wird also zweckmaBigerweise folgenden Rechnungsgang einhalten: man lést zunachst 
Gleichung (25) fiir e’ mit den dem speziellen Problem entsprechenden Rand- und Anfangshbedin- 
gungen. Diese bestimmen im allgemeinen ¢’ nur bis auf eine noch zu ermittelnde Funktion. Darauf- 
hin lést man Gleichung (24) fiir h’ mit den zugehérigen Anfangs- und Randbedingungen, wobei 
eine weitere noch unbekannte Funktion auftritt. Mit den gewonnenen Darstellungen fiir e’ und 
h’ geht man nun in Gleichung (30) ein und erhalt Darstellungen fiir die Komponenten von w’. 
Damit gewinnt man gemaf (20) und (7b) auch Darstellungen fir alle Verzerrungs- und Spannungs- 
groBen. Die damit noch zu erfiillenden Randbedingungen und (22) reichen hin, um die noch offen 
gebliebenen Funktionen zu ermitteln. 


5. Anwendungen. 1. Einachsiger Verzerrungszustand. Ein einachsiger Verzerrungs- 
zustand liegt vor, wenn ein Grundwassertrager, der von unverschieblichen, undurchlassigen verti- 
kalen Wanden eingeschlossen und von einer unverschieblichen, gleichmaBig durchlassigen oder 
undurchlassigen horizontalen Schicht nach unten begrenzt ist, an seiner Oberflache durch gleich- 
maBigen Druck belastet wird. Hierbei mu8 an den vertikalen Begrenzungswanden Schubspannungs- 
freiheit vorausgesetzt werden. Es soll der Fall eines einachsigen Verzerrungszustandes bei un- 
durchlassiger, unverschieblicher unterer horizontaler Begrenzung (Abb. 3) behandelt werden. 

Zur Festlegung der Anfangsbedingungen wird vorausgesetzt, daB die Last P zur Zeit t = 0 
vermittels eines durchlassigen, starren Stempels plétzlich, aber stoBfrei aufgebacht wird. Un- 
mittelbar nach der Lastaufbringung sind allerdings die Beschleunigungsterme nicht vernach- 
lassigbar, jedoch wird sich laut Voraussetzung 4. bereits nach sehr kurzer 

Zeit ein quasistationarer Bewegungszustand einstellen, so daf diese ,,An- 

-laufzeit“‘ vernachlassigbar ist. Da es sich um einen Konsolidierungsvor- 
gang handelt, waren die Materialkonstanten mit dem Indexc zu ver- 
sehen. Er wird jedoch im folgenden weggelassen. 

Fiir den einachsigen Verzerrungszustand lautet (25) 


87e! de! de’ 
ae 48 Gy = Ba 
und (24) 
e?h' ees: 
re ee Os (32) 
Aus (30) wird 
Abb, 3. Schema fiir den ein- orc? 2 qd i) 1 de' oh’ 
achsigen Verzerrungszustand. az Pa rae : Giesae — #18 e — ky % oz? (33) 


worin ¢’(z, t) die Verschiebung des Festkérpers in der z-Richtung bedeutet. Gleichung (22) nimmt’ 
die Form an 


ag! : 
ae — Cale : (34, 
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Zur Lésung von (31) fihren wir die Laplace-Transformierte 
co 
BR Mel}. = [Sigel dt (35) 
0 


ein. Zur Zeit t = 0 + ist nach Voraussetzung 4. die quasistationdre Grundwasserbewegung bereits 


vorhanden, die Porenvolumina fiir z> 0 haben sich jedoch noch nicht geandert. Es gilt dem- 
nach die Anfangsbedingung 


fe Og oe Ore == 0's, (35a) 
Daher ist 
. od! Paes (36) 
at ; 
Somit lautet (31) im Unterraum 
EE diy = 
Ga 6G, *sk=0. (37) 


Die zugehGrige charakteristische Gleichung 


w—xngu—xus =0 (38) 
hat die Wurzeln 
Myre =P (38a) 
mit 
ais) 
p= % (38b) 
und = oo 
= V (8) +s. (38c) 
Die allgemeine Lésung von (37) ist daher 
E = C, 6+ 24+ €, &b—o)=, (39) 
Mit 
ee) 
He UN eel hes dt (40) 
0 


ergibt sich aus (32) bei Verwendung von (36) und (39) 


fH + w)z —W)s 
ae =, (a = 4 Cy eB—%) 4] , 
Daraus folgt 
Be Cees G@  g—a) 1 Cig 4] 
—~ k 1@+ op © Ipsos an 7 Miata 4 ( ) 


Es gelten die Randbedingungen 


dH (42) 


Damit wird aus (41) 


H=% Fpl tgp —u—|, Sx det 0! tga al (43) 
Mit * 
Z= RC} fo evds (44) 
ergibt sich aus (33) 
aa Bee <r Hin ae —2pE—kyn 
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Daraus erhalt man mit (39), (43), (38), (38a, b, c) 


dZ Sits ees 16 i) C, I Cea 
% _ Ceo +o)e 4 C, 0-7 Ts pear. Gea z 


1—2y + w 
C, G C 
GL oy (Go oe? s 
Aus (34) folgt andererseits mit (39) 
= = f = C, e(B + @) 2 |. OF e(B—) 2 : 
Zz 


Durch Vergleich von (45) und (46) erhalt man 


G_ b+o)! Cae (to) ey 
aT apy ae 


und 
C, 
(8 + w)? 


Cy 
(6 — wo)? 


iC, 20s 


Aus (46) folgt durch Integration 


C r C, Seay 
A= eet ae ae Je + C,. 


Wegen der Randbedingung zs = 1... Z = 0 und (47) erhalt man schlieBlich 


yee aa ef +0! Sin [w (I—2)] . 
Aus (10) folgen wegen €,, = yy = 0 und ¢,, =e’ die Beziehungen 
i} a , =e yp , 
OS Bd = te (le whl 2-4) © 
und 
eS l—yp - 7 
eT ET ETS 
Ist 


iz 
“Gees 


die Belastung des Stempels je Flicheneinheit, dann gilt wegen (51) 


1—vy 


1= a0 tna =I &)=9- 


GAabee <a e 
Gleichung (53) lautet im Unterraum 
E.-6 i OAL <9) = 29) g 


1—y nee 


Damit ergibt sich aus (46) und (47) 


we q “4(1+ 7) (1—27) (B+ @) ea! 
: 2s (1—») (Gof ol —B Sin wl) 


und damit aus (49) 


74 Xt) (1-2) ef Sin [ol —2)] 
s l—y oCoj ol—B Ginwl~ 


Die Riicktransformation erhalt bei Beachtung von (38b) und (38c) die Form 


e+10 
¢’(z, t) _ 9% (1+ ) (1 — 29) oie at ; | est Sin [VP + x s(l— 2] We 
eae ami Js Vp + x sCoj (Vp? + xsl) —B Gin (VP + ws) 


Durch Einfiihren der dimensionslosen GréRen 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 
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und der Substitution 6? + %s = s, erhalt man nach Weglassung des Index 1 


c’+i0 
= 1+) (1—2>) ee é™ Gin[Vs (1—2*)] 
em ef lat—preee 1 ee ee 
Soret 3 2m (Vs Cof Vs — B1 Gin Ys) (s — B22) : (8) 
c’—1ic 


Man erkennt aus (58), daf die Abhangigkeit von ys nur scheinbar ist. Es kann daher die Lésung 
mittels des Residuensatzes gefunden werden. Die Nullstellen des Nenners sind gegeben durch 


S27 (59) 
und durch die transzendente Gleichung 
= ‘8 
Tqyjs = BL" (60) 
Gleichung (60) besitzt nur dann fiir /s eine von Null verschiedene reelle Lésung, wenn 
1 7B 
Bi < py: ? 
d. h. nach (38b), wenn 
tgl>2y2. (61) 


Zur Abschatzung der Gréfenordnung von x,g gehen wir von (31) aus und betrachten den 
Endzustand (t = oo) bei Vernachlassigung des Auftriebsgliedes 


de’ 
48 Fee 


Da nunmehr in jedem Querschnitt (0:,);— « = — q = konst. ist, ist (e’),— . im ganzen Bereich 
konstant. Es gilt daher nach (34) 


Cisco = AE )raw $C. 
Weren -z.=1...(C ),- ~ = 0 folet 
| (2 = 0 =—(l—2) (6) mo 
und an der Stelle z = 0 
(0): © = S41 ee 


Zz 


I tl 


Daher ist mit (53) 


0" =— (Cinco = PEE 9 = — my) 09 (62) 


Hierin bedeutet a,, den von Terzaghi! eingefiihrten Verdichtungsbeiwert. Es ist also 


meee Y) (Ll — 2.7) 
“oc “@ —v)(1—™) * Ne 


Demnach folgt aus (26) und (63) 


1 § = Ay, (1 — m9)? (7K — Ym) - | (64) 
Fiir blauen sandigen Schluff-Ton (Wien) erhalt man die Werte? 
ty = 0,6. a,, = 3,0 - 10-? cm?/kp , 


Ve = 2,1*10%kp/em?, yy =10-*kp/cm’. 
Damit gibt (63) 
ye 1,84107> cm. 
Die Ungleichung (31) ist in diesem Fall nur fiir 1> 1,54 km erfillt. Wenn, wie wir voraussetzen 
wollen, die Dicke der Tonschicht geringer ist als 2 /2/x,g, besitzt Gleichung (60) fiir ys keine 


1 Siehe S. 270 des Buches von Terzaghi-Jelinek, Fu8note 1 von Seite 225. 
2 Gnosgritie Wasserkraft u. Wasserwirtschaft 33 (1938) S.5, speziell die Tabelle von S. 9. 
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von Null verschiedene reelle Lésung; hingegen existieren unendlich viele imaginare Lésungen. 
Wir setzen 


/s; =i9;, (65) 
worin 0; positive Zahlen bedeuten. Aus (60) und (65) folgt 
tg 0; = 5 (ee oe (66) 
Die zugehérigen Werte s; sind aus (65) 
5, = —-0, Gf = "0, 1203's ae (67) 


Mit Hilfe des Residuensatzes findet man schlieBlich aus (58) nach einfachen Umformungen 


a €@ 


2 22) 4* 
pe ets dnan Arlee rie ay See sie (68) 


ae j- tos... (+ BP) (of —B 1+ BP) sino; 
Ist 
pl= ae cae (69) 
dann erhalt man aus (68) durch den Grenziibergang fb— 0 den Ausdruck 
ce) ; —j+1 =e 
q% (1+ ») (1 — 2») 8 (ie late : pene 
ct = ees {a 2*) 72 Ze Sy eee, sin i go re Gs Fo H (70) 
J=U,1,2... 


und daraus die Setzung an der Oberflache (z* = 0) - 


—__.,— |. 71 
1—y Le me (27 + 1) ( ) 


(C*),« pen Oe?) ; 8 Ss ae 


Diese Gleichung stimmt mit Terzaghi-Fréhlich! iiberein. 

Fordert man wegen der Ungleichung (69) z. B., daB 61 < 0,05 sein soll, so erhalt man fiir 
die obigen Zahlenwerte die Forderung 1 << 55 m. Auf diese Weise kann man den Fehler abschatzen, 
den man durch Vernachlassigung des Auftriebsgliedes begeht. 

Es soll nun noch gezeigt werden, daB sich beim einachsigen Verzerrungszustand auch fir 
0,, eine Differentialgleichung finden la8t und daB sich die Differentialgleichung (29) fiir h’ durch 
eine zweiter Ordnung ersetzen ]aBt. Aus (31) folgt mit (51) 


ee 
Of Ae ape = oe (72) 
Aus (6b) mit (17b) erhalt man zundchst 
a0, op’ , 
<p fe ee ) Vy) © (73) 


Wendet man darauf den Operator 


an, so folgt wegen (31) und (51) 


2 / ep’ 
8 (0 meni 0 es 
bzw. wegen (5b) auch 
@ / eh’ oh' oh’ 
alas "855 — pr) = 9 (75) 
oder 
2h! t ' 
Fae ee =F, (76) 


1 Siehe FuBnote 2 von Seite 225. 


XXX. Band 1961 G, Heinrich u. K. Desoyer: Theorie dreidimensionaler SetzungsvorgangeinTonschichten 235 


worin f(t) eine unbekannte Funktion der Zeit ist. Man uberlegt sich leicht, da8B f(t) auch fiir den 
Fall einer zeitlich konstanten Last tatsichlich eine Funktion der Zeit sein muB: fiir eine von der 


‘Zeit abhangige Belastung q(t) je Flacheneinheit gilt laut (35) 


i 1—y , 
(6525 26s Eien) 35) (e’)z-0 = — q(t) fir t>0. (77) 


Demnach ist wegen (32) 


| oS OE 8) = 2): dq(t) 
dz }, 9 ke(i—»)  - dt * 


Wegen der Randbedingung (h’),_) = 0 fiir t> 0 gilt a = 0 fiir +> 0. Demnach folgt 
aus (76) O Jz =o 


=f(t). (78) 


a (1+ v) (1 — 2) dq(t) oh! 
kp (1 —) "weet (ze 


Mit Beniitzung von (4b) ergibt sich daraus 


Aeeeo a (1 + v) (1 — 2») dg(t) 
F(t) = 3 (s)e=0 Raia) ae (79) 


Aus Gleichung (79) erkennt man, daf f(t) im allgemeinen tatsichlich von der Zeit abhangt; f(t) 
verschwindet dann, wenn q zeitlich konstant ist und wenn tiberdies auf eine Beriicksichtigung des 
Termes mit x, (Auftriebsglied) verzichtet wird. In diesem speziellen Fall gilt 


e2h' oh! 
age ee ct (80a) 
oder wegen (5b) 
ap’ dp’ 
ee ag (80b) 


Dies entspricht der in den Lehrbiichern tiber Bodenmechanik bisher ausschlieBlich verwendeten 
Ausgangsgleichung fiir die Berechnung von Setzungsvorgangen. 
Verschiedene Autoren’ versuchen, diese Gleichung fiir den raumlichen Fall in der Form 


Ah’ —x — =0 bzw. Ap’ —x 1 0 (81a, b) 


zu verwenden. Wie (29) zeigt, ist dies auch bei Vernachlassigung des Auftriebstermes unzuliassig. 
DaB (81a) und (81b) auch fir praktisch wichtige Spezialfalle nicht erfillt sind, kann leicht an 
Hand von Gleichung (160) festgestellt werden. 


2.Setzungen von Tonschichten groBer Machtigkeit bei zylindersymmetrischer 
Lastverteilung. Ein Grundwassertrager, der als unendlicher Halbraum angesehen wird, werde 
an seiner horizontalen, iiberall durchlassigen Oberflache zur aint) 


Zeit t = 0 beginnend, durch eine stoBfrei aufgebrachte zy- 
lindersymmetrische Druckverteilung g(r, t) belastet (Abb. 4). _—<ciiin, 2 
Dabei wird 0q(r, t)/Ot = 0 vorausgesetzt. Es soll der zeitliche hee ag 


Verlauf des durch die Belastung eingeleiteten Konsolidierungs- | 
vorganges ermittelt werden. Wie im vorigen Abschnitt ge- 1 


zeigt, spielen die von der Schwerkraft herriithrenden Terme Abbe iiacanelcische Belactune 
des Halbraumes. 


in den Gleichungen (25), (28), (29) und (30) schon fiir den 


- einachsigen Fall nur eine untergeordnete Rolle. Fiir den 


zylindersymmetrischen Fall, bei dem die StérgréBen mit der Entfernung von der Belastungs- 
stelle abklingen, sollen diese Terme zunachst vernachlassigt werden. 
Die Gleichungen (25) und (30) lauten dann 
fe al 
Ne (82) 
und 


; 2(1—») fo 
De a ey Oc a (83) 


1 Siehe FuBnote 3 von Seite 225. 
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mit 
ence ee (24 
und : &s % 
or = = - > Z ae i ee 2 (sa) 


worin ¢,, Cy, ¢, die drei Einheitsvektoren in den Richtungen r, @, z (Zylinderkoordinaten) bedeuten. 
Es gelten folgende Anfangs- und Randbedingungen: 

lstire te. pgson Ve pigtet == ie (86) 

Dhar. ea p=) Seah’ 205 (87) 

a the oy a0) 2=07..¢,=0, (88) 

Ay furs = 0, z= 0 x. 0;, = — q(r.t) « (89) 


AuBerdem muB gefordert werden, daB alle StérgréBen fiir z+ oo verschwinden. 
Zur Lésung von (82) fiihren wir 


co 
I os Qe =f ened (90) 
0 
als Laplace-Transformierte von ¢’ ein. Wegen der Bedingung (86) folgt damit aus (82) 
Ak =o s.E. (91) 


Zur Lésung von (91) benutzen wir die Hankel-Transformation. Fiir eine Funktion f(r, z) lautet 
die Hankel-Transformierte (siehe z. B. 1) 


ice) 
9, {f(r 2)} = J F(t. 2) 7 J,(r p) dr, (92) 
0 
wobei J, die Besselfunktion y-ter Ordnung bedeutet. Die Riicktransformation lautet 


f(r, 2) = (f(r. 2)} p Ir p) dp - (93) 


AuBerdem gelten die Beziehungen 


bo (Afr. 2} = (Za —P*) Soll 2)} (94) 

und 
: (4 — a] 0. 2} = (ia —P*) Gil 2) (95) 

Mit 


co. 
Ho{E} = f Er Jy(r p) dp 
0 
liefert (91) bei Verwendung von (94) 
d? : 
(jz —P° #8) So(B} =0. (96) 
Die allgemeine Lésung von (96), die fiir z—> 00 verschwindet, ist 
Ho {E} == A(p, s) e— Vp? t+usz : 
worin A(p, s) eine willkiirliche Funktion bedeutet. Die Riicktransformation nach (93) liefert 
ee) - 
B= f A(p,s) VP" p Jor p) dp es 


* C. J. Tranter, Integral Transforms in Mathematical Physics, S.2 und 12, New York 1959. 
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Die Elimination von @e‘/@t aus den Gleichungen (24) und (82) liefert mit (27) 
glee PE 29) 


is Oy eAna == As’. 
Setzt man 
1 me 
Q a ( a (1 2 yy h’ =s TI ; 
—y 
so folgt aus (98) 
A(H — E) =0. 


Nach (94) ergibt sich mit 
(oe) 
Do{H} SER Jo(r p) dr 
0 
aus (100) 
d2 
fei ) ){H — E} =0. 


Somit ist 
Do {H — E} = B(p, s).e—P* 


bzw. nach der Riicktransformation gem&8 (93) bei Verwendung von (97) 
© ioe) 
H= f A(p,s) eV ****p Jor p) dp + f B(p,s) eP* p Jor p) dp 


Hierin ist B(p, s) eine weitere willkiirliche Funktion. Aus der Bedingung (87) folgt 


A(p,s) + B(p, s) = 0 
und damit 


H= f Alp, s) [e—Ve*+%s*_ 9—p2| P Jo(r p) dp 
0 
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(98) 


(99) 


(100) 


(101) 


(102) 


Bedeuten u bzw. w die Verschiebungen des Festk6rpers in radialer bzw. vertikaler Richtung 


uo=u Cp we 
dann liefert (85) bei Beachtung von (27) 


oe 1 
A ee ap +20 +) dare & Oz 
und 
ge? it 62 
ee agahe= —.¢ +20 (1+) dove 5 - 
Mit W = L{w} und U = L{u} 


wird aus (104) und (105) bei Benutzung von (90) und (99) 
Ton eb 2 9) oe 


AW 7 3 0s 1—2v dz 

und 
1 lL O6E , 2(1—») eH 
eee ay ot as ér 


Durch Verwendung von (97) und (102) sowie der Beziehung 


2 [Jolt P)] = —P Slr P) 
erhalt man aus (107) und (108) 


AW = [40.9 POP eer — |p pase VP FM p Jr) dp 
0 


\ 


(103) 


(104) 


(105) 


(106) 


(107) 


(108) 


(109) 


(110) 


938  G. Heinrich u. K. Desoyer: Theorie dreidimensionaler Setzungsvorgange in Tonschichten _Ingenieur-Archiv 


und 


co ‘ 
(4— ss) U= | p Atos) pags Yee He) ae (111) 


Wendet man auf (110) die Hankel-Transformation §), an, so ergibt sich wegen (94) und (93) 


2(1 —») 


(a *) S07} = A(p,s) = pe—P*— pt + ns ele txen) ay) 


Analog ergibt sich durch Anwendung der Hankel-Transformation , auf (111) mit (93) und (95) 


a? © 2(1—%) ps __ 9 -Jp? use 
BS *) {0} =P Alp.s) pa, e—P+ —e—VP?4 |. (113) 
Die allgemeinen, fiir z—> co verschwindenden Lésungen von (112) und (113) lauten 
ie—s 1 3 
9 {W} =— A(p,s) hi = gee 4) ph fuse WPT] 4 C(p,syeP* (114) 
und 
is Ls 0 pipe PS eee: : eet 115 
rely Ae 3) Oo ee ee + D(p, s)e~?*. (115) 


Hierin sind C(p,s) und D(p, s) zwei weitere willkirliche Funktionen. Die Riicktransformationen 
nach (93) liefern 


co 


ae ; 1 
Ww —— | atp.s) | ayy cer +g (pees VPs) — C(p,s) ol pIuirp) dp (116) 
; 
und 
ice) 
= Des” — pz Po o—\p?+usz —pz 
ee a, Fo OE — D(p, s) e~P*, p Jy(r p) dp . (117) 
0 
Gleichung (22) gibt bei Beachtung von (85a) und (104) 
a) 1 ow b 
Geeebe (118) 
bzw. im Unterraum 
é i ow . 
Gt >) U+SE=E. (119) 


Setzt man in (119) die Ausdriicke (117), (116) und (97) ein, dann erhalt man bei Verwendung der 
Beziehung 
i) 1 
e ae = J,(r p) = p Jr p) 
die Bedingung 


(ee) 


| p2@9)— Fay 40.9—P co. 8)] Pp Jolrp) dp = 0. (120) 
| 


Da diese fiir alle z > 0 und r = 0 erfiillt sein mu, folgt daraus 
esi . 
p D(p, s) — ah eae A(p, s) —p C(p, s) =0. (121) 
Zur Erfiillung der Randbedingungen (88) und (89) driicken wir zunichst die Spannungen des 


Festkérpers durch die Verschiebungen u und w aus. Wir setzten (21) in (10) ein und erhalten 


1 


/ 1 7 / RS f 
Sea =a (Vou +0, 7) +555 eG]. (122) 
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Mit (85) 1aBt sich (122) bei Beachtung der Zylindersymmetrie in der Form schreiben 


, , , ou v 1 /éw du 
(oy to} o i) ae ——— ¢' = 
(oi, ae, ae Gace Poe 0 aes 
7 , , i! u Vv ‘ 
Ce eC C9) (ee 0 = Sond 
gr ep Vp o (1+ ») r a spear rae 
; , ; Low Ow ow v 
ERD Raa eae |e ee Seer gl 
2 ¢, = | ates 0 ee ey © 
Die Randbedingung (88) lautet gemaB (123) 
é au 
| 5 Sik 
or 62} ,~9 
bzw. im Unterraum 
oW aU 
ce a ae : 


Unter Verwendung von (116), (117) und (109) ergibt sich aus (125) 
ie) 
| 


Da (126) fiir alle r > 0 gelten muB, folgt 


2 Sfp ere 1 — 
(co VP Fs — 7) A. 5) —p Cp. 8) —p Dip. s)] p Slr) dp = 0. 


1—yvy 


2p 
(Bo Pt tes Tay) AP. 8) —P Cp, s) —p Dips) = 0. 


Die Randbedingung (89) driickt sich mit (123) in der Form 


1 ow v r ez, 
Tee Te 
aus, bzw. 
1 OW v 
sce) ee brag El, a 
mit 


q Or, 8) = L{glr,2)} Sih o—" q(r,t) dt. 


Mit Verwendung von (116) und (97) folgt aus (129) 
ie.e) 


[| 4@-9—p 6. 9)] p Ilr) de =—a 01 +9) Or9). 


0 


Es ist zweckmaBig, darin Q(r, s) ebenfalls in Integralform darzustellen. Dazu setzen wir 
es) == 
Do {Q(r; s)} ui Q(r, s) r Jo(r p) dr = Q(p, s) - 


Dann liefert die Riicktransformationsgleichung (93) 


lr. 8) =f Op.) p Julep) a 
Aus (131) und (133) ergibt sich 


1 ee A(p, s) — p C(p, s) +a (1 + ) Q(p,s)| p Jo(r p) dp =0. 


0 


Da diese Beziehung fiir alle r > 0 erfiillt sein muB, folgt 


7 fe A(p, s) —p C(p, s) = —« (1 +) Q(p, 8) - 


SE a ee we 


(123) 


(124) 


(125) 


(126) 


(127) 


(128) 


(129) 


» (130) 


(131) 


(132) 


(133) 


(134) 
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Die Beziehungen (121), (127) und (134) bilden ein inhomogenes Gleichungssystem fur die drei 
Funktionen A(p, s), C(p, s) und D(p, s). Sie liefern 


a(ltv)(l—2»)u%s Ole (135) 
Alp, *) = aha Q(p; s) 
Che _a(1+») [(1—2») p VP pet xs—(1—»)xs] QP. 8) (136) 
P: (1 —20)p (p—Vp?+xs) + (1—v)xs iB 
i jpetus = 
D(p, 3) = ee eee (137) 


(1 —2») p(p— Jp? xs) + (1—»)xs 
Durch Einsetzen dieser Ausdriicke in (116), (117), (97) und (103) erhalt man 


/p?+usz 

[a —»)xs(1+ pz) —(1—2») p |p? a P* 4+ (1—20)p Vp?+ xs +e se Vp 
W=o(1+ 2 sh ees 
ce »| (1—2)p(p—Vp? + #8) + 1 —n)xs 


x Q(p, s) Jo(r Pp) dp (138) 


—|p?tusz _ 
19 f Wane R* ja(e oy ip ems O(p. s) p Jy(r p) dp ; 
(1 —29)p(p— |p? tus) + —»)xs 
(139) 


xs eo Vpt+xs % = 
Q(p. s) p Jo(r p) ap ; (140) 


(1 —2) p(p—p? + %s)+ (l—»)zs 


E=—oa(l sci 


pe allt eo 
H =—a(1+ »)(1—2») | Ge aa Q(p, s) p Jo(r p) ap - (141) 


(1—2) p(p Vp? +> + (1 


Die Gleichungen (138) bis (141) stellen die Lésung des Problems im Laplaceschen Unterraum dar. 
Die Riicktransformation kann erst durchgefiihrt werden, wenn man sich eine spezielle rotations- 
symmetrische Belastungsverteilung q(r, t) an der Oberflache vorgibt. 


6. Lésung fiir zeitlich konstante Gleichlast. Im folgenden sei 


jee gq =konst.c fitr—aa, (142) 
0 finn se 
Damit erhalt man aus (130) 
q ze 
as f f 
66a) urr<a, (143) 
0 fury >'a, 
und aus (132) folgt 
Qlp.8) = | fr Jolrp) dr = 22 Jap). (144) 


= 0 


Dieser Ausdruck ist in die Gleichungen (138) bis (141) einzusetzen. Wir formen die entstehenden 

Ausdriicke durch Rationalmachen des Nenners und darauffolgende Partialbruchze1legungen um 

und erhalten: 
ee) 


W=« (1 +)qax| {2 


0 


Sa (» = Vp? Fs) (eps — el + #82) 


(1—y) (1 +pa(l a 


+ = 2eP 
x? s2 


(p + Vp? + s) (e-P#— eV? +82) 


ci ye oe 
Se —2y) p? a al - */p? 4 x5) a + pz—v)e—P* 4 ye—Vp tHe 


ee Jo(r p) dp : 
(145) 


VON 


KX a \ ’ ‘ r ‘ ‘ + ‘ i 
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1 [Ql —»r)z 
Uso (1 + 9)gax | [ ( =) (p 4 pt pws) o-P* 
0 

a 1 F Il—y (SD ( —pz 1 Menon 

a e \p t xs ve + ( —v)e Fam) 

i= = a 

ae (P+ Vp? + 8) (oP — eV Fe) 

; 1 Loa ? 
geG “Fj (ke : Lp »—( tr) lp? +xs +x] 
x [(ps—9) PL 9) IP I (ap) Julep) dp, (146) 


B=—a(l-+»)0—29)¢0x [2 (1 4 Ee 


1—y —_ —— 
1 Cae) pt (Ly AG P ; \P* xs} e—Vpt+xs= Ji(a p) Jo(r p) dp , (147) 


ee) 


So (1 +0 —20)90% | (p + Vp? +xs + x8) 


i a r 
aye ee (1—») |/p? pet xs]}(e- p2__ @—Vp? +xs Free) MEP) JOOP) gy (148) 


Fir die Riicktransformation benétigt man 


—Vpi+xsz 1 ae eae ee 
gals ce — oP Exfe ( ma | = L¢ —P Exfe()/* Tees = ; 39) 
—|p* + x82 l PE RT ee ae 
ee Ais t Zz . nz pit t Zz $ prt 
: ae =e (st xploreeel/ e+) +(e) eree VV), 
(150) 


Q-1 Sarge iz 1 eal eres f prt 
(L—2») p?+ @—»j'xs)  2¢Q—vp° : aaa Pa | 


—Vpi+usz 
\p? +xse- 2 P 
pfeil 2 


“Ss 
Sy (ra ale 
a Gate) (152) 
-yJVP? + %s er eee = Tau iad pee e) 5 | 2? yes) 
go pele 2” p Fae cool ky eP* Erfc le ae 
fet yt 
a ak Bate (| ule Va c ), (153) 
\pr ene BS = es era 
. | Ode Caged Vind THe | Ia a *) 
ce as = re oe e rfe\~ 4¢ '1l—vV x 
v ss aso wz? pit 
ne aes /=) : (aE) 154 
aoe Efe ( At l1—vy x # (1 — vy Jae 7 >. (154) 
worin a | 
Erfo(x) = a als e-* dv (155) 
‘ A 


die normierte komplementare Fehlerfunktion bedeutet. 
17 
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Die restlichen erforderlichen Ausdriicke erhalt man durch Spezialisierung, wenn man in (149) 
bis (154) z = 0 setzt unter Beachtung der Beziehungen 
Erfc(0) =1, Erfe(oo) =0, Erfe(x) + Erfe(—x) = 2, Erfe(— Co) == 2/5 
Damit erhalt man aus (145) bis (148) bei Beachtung von (90), (99) und (106) 


== (1 +oga [{em[a—ay® Pinte = + Erfe (“2 aes 2) +2(1—») pz 
0 
+2 (8—4 9 + 49%) — + (2 (1—2») ps—(1—4 »%)) Bovfe (\/2* = 
1—2y ptt ‘ 
ae Moe fs pit ia (l—»)?_ x v pit 
—(a ”) pets a(t 2” +-292)) Bfo]/* +yv(1+pz re. * Ble (2, *t) 
1 —2v1/p?t ee 1 pz a ad p?t 
=. iE ak, jails, Efe (|/** = 
y 1—2»p pt 
i=;t GF se m2 Y pit 
+e Berto (|/24 re ae 5 


(“+ 3 Poids t 


olay “| eon Lien) Ti ap (157) 


(156) 


1—2y re 


\x oe 
ua (l +9)g0| fer] —25) (ate YE 4 aan F\_2)+20—»)ps 


—La—4y—La—an eps +1—2y Bete (Ve) 


1—2y prt 


—(a—» ps—h. a —29%) rinto |/t +» (pz—)e_ G— * Hig. te Zaye 


2 f 
1—2y 1/p?t —P- 1 ( 2 j= 
= _— ——\@ x i eP* Ney ae 
jz = |+ re Erfc re _ - 


+ 


Pall ced? v bs 3 —- 
<n (ls) FP ets (VEE yt) 
——_ 4 
1 22 ba 
1—2 1 /ptt m —(s+e) J,(@ p) Jy(r p) 
Le eee pe ee ; me 


1—2» pt y 
6 Gaye iS 


+. 


l—vp t 1—y 


ee) 
ee 2S eer Brfe =) 
x 2 
6 


ae 1—2y pit ——, 
tye (—) x Exfo( : y= 
1—v % 


1—2p p*t v 


Be var Seen er we pt 
| ve ( ) 1 Efe ( er + Ve h J,(a p) Jo(r p) dp * (160) 


I— 


erte (V2 + y=)| Jap) Jory) dp, (159) 


Man iiberzeugt sich leicht, daB h’ die Bezich . a : 4 
Pigeide Dowereatey e Beziehung (81a) nicht erfiillt [vg]. die nach Gleichung (81a) 
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Aus (157) und (158) ergibt sich, da8 bereits unmittelbar nach der Lastaufbringung (t = 0-1) 
von Null verschiedene Verschiebungen w und uvorhanden sind im Gegensatz zu dem vorher be- 
handelten einachsigen Problem. Es ist 


(wo = ol +9) Ga f (I <2) eps AEE) IP) ay, (161) 
0 

(uo =o (1 +9) ga | ze-P* J,(ap) Jy(r p) dp . (162) 
0 


Fiir die praktisch besonders wichtigen Werte der Verschiebungen an der Oberflache (z = 0) erhailt 
man aus (157) und (158) 


ee) 


ite}, 6.0 (1 +n 0—nae | [2—ner oe 


0 


1—2» p*t eee 
Vig eh (=a) Sra v at id AG Ji 
rere oa Exfe (-7—/® ] eS ac! 
co 
(,-9=—a(+y)40| |a—2)—a—») Bere y= 
t>0 6 as 
__ 1—2» pit 
+ ve oF ite (y=) ce peney dp . (164) 


Aus (161) und (163) bzw. aus (162) und (164) erhalt man tibereinstimmend die Anfangswerte an der 
Oberflache 


ioe) 


lim (w);—0 + — lim (w),—0 W(X) (1 + 9) qa | Ji(a p) Jol" P) ay, 
z—0 t>0+ F P 
zal +9) go (2) fir r<a, 
< 2 (165) 
a (1+ y)q [e(+)—(t —%) (+) fir ra, 
lim (u),-04 = lim (u),-0 = 0, (166) 
z—>0 t>0+ 
worin 
? $$$ $$$ 
E(é) 4, /1— & sin? y dy (167) 
- und 
2 d 
sbi abht SN 168 
KO = | Sears (68) 
0 
_ bedeuten. In gleicher Weise erhalt man aus (159) und (160) 
(e’).-04 = 9; (169) 
Z2>0 
(170) 


CO 
(Joy hooe = gas o?* Sy(ap) Jo(r p) a 
; z>0 
r 17* 
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und 
(ce) 
(e’)s-0 =—x (1 + ») (lL—2») q of [pee 
t>0 
1—2» p?t 3 
) ~ (L—»)? «- t 
+5, pe ake Oe Fen J,(r p) dp, (171) 
Parpees 0 (172) 
t>0 
Aus (169) erhalt man 
lim(e' 6 4 == 0. (173) 
2—>0 
jedoch aus (171) 
= a(1-+»)(1—29) | ¢. 
1 eens — = giurr<a, 
lim (e"),-0 =F PC ga | ap) Solr) dp = ae (174) 
oe 0 0 fir r>a. 


Aus (170) bzw. (172) erkennt man 


ee) 
“ : fiir™ ri @, 
lim (A) Viv’), =0 4 = go | J\(a p) Jo(r p) dp = i : (175) 
x0 ‘ REL ree txs 
hingegen 
lim (Ay yw h’),=0 = 0. (176) 
t>0+ 


Man sieht hieraus, daB ¢’ und h’ zur Zeit t = 0+ fiir z = 0 eine Unstetigkeit aufweisen. (169) zeigt, 
daB die Anfangsbedingung (86) erfullt ist, aus (174) erkennt man, daB zur Zeit t = 0+ an der 
Oberflache fiir r < a ein ¢’ vorhanden ist, wie es einem einachsigen Deformationszustand entspricht 
(vgl. Gleichung (79)). Die Gleichungen (172) und (176) zeigen, daB® die Randbedingung (87) erfiillt 
ist; aus (175) erkennt man, da fiir t = 0+ unmittelbar unter der Oberflache A, yyw h’ gleich der 
Belastung q ist, wahrend A) yyw h’ laut (170) mit wachsendem z gegen Null geht. 


Bei Beachtung von 


lim ue 2 zm Bixfe( ae Ea ) Erte a I 25 (177) 
erhalt man aus (157) bis (160) fiir den Endzustand (t = oo) 
z 
(waco =a +9)gaf [pe +2 —s)] PAE ICD) gp, (178) 
; lee) 
(na =a(l-+nga{ [p2—G—2)] cP MEP TED ay, (179) 
te eee 
Pe a(t (1 2.5) qa f P J,(ap) Jo(r p) dp, (180) | 
ov Bh’): 0 = 0. (181) 


Die Gleichungen (178), (179) und (180) gehen, wenn man den Grenziibergang zu einer Einzellast P 
durchfiihrt, in die bekannten Gleichungen von Boussinesq iiber. Dazu setzt man’ 


P 


tS (182) 
und fihrt den Grenziibergang a— 0 bei Beachtung von 
4, a2 P) P 

care a cae “2. (J 83) 


durch. 


poa\ ae 


ey eee ye 
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Aus (178) und (179) erhalt man fiir die Verschiebungen an der Oberflache (z= 0): 


(w)s = co =2a(1—*) ga | AEP) ACP) dp 
2z2=0 ‘ 
zo (lL) qaE (2) firr <a, 
as ; (184) 
(1 — ») qr B(7)—(1—3) (4) fiirr=a, 


—=(1+»)(L—2») qr firr = 4.5 


(u)inco =—a (1 +9) C=2)) 0a) Eee dp = 
z=0 ; I — (1+) (1—29)q> firr=a. 
(185) 


Die Spannungen im Festkérper ergeben sich aus (123) mit (157), (158) und (159). Man erhalt bei 
Beachtung der Beziehung 


ae Wale) =P |Jolr p) — 0 (186) 


die Ausdriicke 


vin= aa | fe ra 29) (ee | ae #\_2)420—mpe 
0 
1 1 : Nee [pre 
a he (Ee (2 pz 1 — 2) Exfo( “ee 


1—2y pt 


1 2 t == (V9) oe v prt 

—(@—»)pz— 0 —2»)) Erte VE + »(pz—n) Oat Brfe(72 
Pore ES 2 

1—2y]/p*t — % ieee ne pt 

ef yee [+ ternete( aS 

29 pik- te ie 

ptt 

Bee (|/* eee) =) 


ae 
ie) 
Z — pz ne pee 
Jo(r p) — ee dp — cote P Efe ( Tre F*) 
0 


ps 


_ (et 
ae ee pt e (3 4 I rap) 


ia) 
ae als v KA = ] 
2 it x I» % 2 y p? 
cS _—e (a=) Efe (|/%* 2?) J,(a p) Jo(r p) ap , (187) 
Op = 0, - (188) 
~ Shed 
2 2 
cin i | fore] 2(1 —20) 2 ‘(2—Exfe ae eens Er fo |/Z - ‘\—40—»)ps 


+- oe + (1—2») (p2—+ (1 +29)) Bxfo(|/2" 24) 
+ (2 (1 — ») pz — (1 — 20%)) Erfe yet 


tee 2 2(1— 2»), /p?t Et 
+2(y—pa)e ie Exfe(, P tere e 


—y ice \z ae 
qo + pe % 2" v pit 
tees a \— 2 ve 5 ) : Efe ceo 
VBE og 
+ 
ee ee gla Nien I(r p) dp ; sed 
)z % 
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Coy = = 4 = eae mit uw gema® (158) und e’ gemaB (159), (190) 


canna [{erfa— 2») Bt (2 — Eate(/ 5 see Pa) Erfe|/®")—2(1—») ps 
my ee (a —2») pz —4. 6 — 498) Erte (/2* = P) 


1—2y» a 


+(a—»)ps—} a —208) Erte /t —v(pe— no aap * Efe (” are =e 


ptt 
1—29]/p?t ~ x ies8 % 2 prt 
tC oe eee lene Erfo(/4* +E 


1—2y p*t v 


she pz 2 
Ey ot ae hgeaet 4 Ae 
ve Barfo (//<* gag =) 7 


1— 27, /2t — ee 
———-|/F- e Cae I J,(a p) Jo(r p) ap - (191) 
|x x 
Wahrend im einachsigen Fall der Festkérper zur Zeit t = 0+ fiir z > 0 noch keine zusatzlichen 
Spannungen aufweist, erhalt man hier fiir t = 0+ die folgenden Ausdriicke: 


é 
(or denos = 40 | pee? Jap) |lrp) — 9 |dp (192) 
und daraus ‘ 
lim (G7) 04 = 05 (192a) 
(oino4 =—gaf p2e-* Jap) Slr) dp (193) 
on (o7,)h=04 =0, (193a) 
“ 
(Codvnas = 90 { FeP* Klay) A(rp) dp, (194) 
lim (opp)=04 = 0, (194a) 
(6:01 — 9a pze-P* Jap) Jule) dp (195) 
lim (6/,),—04 = 0. (195a) 
Fiir die Oberfliche (z = 0) erhalt man aus (187), (189), (190) und (191) 
, 1—2y pt 
(Orr)=0 =qa | E 1 + Erfe exes cS e UF * Erg, ie =| J,(a p) Jo(r p) dp 


f 1—2y ptt 
=e 2 a rT 
+t2[[a—ao—0—» Erfe|/P! tye OF * Eafe — Ve [seme a, (196) 


q fur r<a,| 


v 
lim (0r))e=0 =— 290 ic Jap) Jorp) dp =} =? (196a) 
0 Os firir Sa; | 
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© 


(Cop)=9 = —a (1 +40 | |a—2n—a—n Erfe 


a 
tye (l=) *Erfe( ye z)| ion i P) g, 
l1—y»y 4 Tp 
1o@) 
—art-+9)qe| |2—Bet ae 
0 
ao et 
v ~ (l=) Bet 
as, * Efe ( ze J,(a p) Jolt p) dp, 
(189) 
: y Be y+ — ly fir sh 6; 
tim (d9»):—0 = —a —— re [ (a p) J(r p) dp = : (198a) 
Oetur srs G; 
~ q firre 
(2)e-0=— 4a f Jx(ap) Jo(rp) dp =} a 
t>0 J . o(r P) dP 0 fir pea ee, 


Die Gleichungen (197) und (199) zeigen, daB die Randbedingungen (88) und (89) erfiilt sind, ferner 
erkennt man aus (196a) und (198a), daB zur Zeit t = 0+ fiir z = 0, r <a ein einachsiger Verzer- 
rungszustand besteht, was bereits auch aus (174) hervorging. 

Fir den Endzustand (t = 00) ergibt sich aus (187) bis (191) folgende Spannungsverteilung: 


(Corawm—va| i (L—pz2)e—?* J,(ap) Jo(7 p) dp + | [pz—(1—2»)] | (200) 
0 0 


EXG):-0 = 0, (201) 


Cia a=—¢ op »f e—P* J,(a p) J,(r p) dp + i [(1 —2 ») — pz] ep Su(¢ P) SiC P) to} eon 
rp 
0 0 


pias = ge if (1 + pz) e-P* J,(ap) Jolr p) dp. (203) 


Wie man sich leicht itiberzeugen kann, gehen die Gleichungen (200) bis (203) fiir den Grenzfall der 
Einzelkraft in die bekannten Formeln von Boussinesq iiber. 


7. Berechnung der Setzungen. a) Setzung unter der konstanten Flachenbelastung. 
Da im Gegensatz zum einachsigen Fall in dem hier behandelten raumlichen Vorgang bereits zur 
Zeit t = 0+ Verschiebungen der Oberflache in vertikaler Richtung vorhanden sind »Anfangs 
setzung“, soll der Ausdruck 


S(r, t) = [w(r, t) — w(r, 0+)], =0 (204) 
als ,,Nachsetzung“ bezeichnet werden. Aus (163) und (165) folgt 


S(r,t) =o (1 +yae| [a—2—a—n Beto 


0 
__1—2» ptt t —— 
+ ye a ee * Efe (7 y=)| Ji(a p) Jo(r P) dp : (205) 


—Y» x" P 


Um einen mittleren Wert der Nachsetzung unter der Flachenlast zu bekommen, soll der Mittelwert 


Sm(t) =p i 2 rae S(r, t) dr (206) 
0 


248  G. Heinrich u. K. Desoyer: Theorie dreidimensionaler SetzungsvorgangeinTonschichten _Ingenieur-Archiv 


als ,,mittlere Nachsetzung“ eingefiihrt werden. Mit 


[ rdlrp) dr == lap) (207) 
0 
erhalt man hierfiir aus (206) und 205) 


co 


Sp(t) = 20 (1 +na{ [a—2) —a—» Beto 


0 


1—2y p*t = P ) 
a eae a V Ana feb 208 
+ve rfe( "= ] 7 dp (208) 
Fiihrt man die Substitution a p = p, durch und verwendet als dimensionslose Zeit 
a? 4, 
so erhalt man aus (208), wenn man statt p, wieder p schreibt, 
ice) 
Sn =20(1 ives] |@ oo Erfe (p //t*) 
Re Lp —\1 72(p) 
i eee ag Erfe(; = -P yi") “2 (210) 
und daraus fiir t = co 
ets 8a 
(S. \ou= Bo vy) 2 9) ne if fe) dp=5"(1+0)(1—-2»)qa. (211) 
0 


Definiert man als ,,dimensionslose mittlere Nachsetzung“ 


eee Sin as 37% 
Sm (Sra)t = 00 a(t) —anyqa>™” 


so wird schlieBlich aus (210) mit (212) 


(212) 


a ee 
eat 3a(1—1) me 9 Gas v = \| Ji) 
rt ea Erfe | — p jt =a. (213) 
6 
pease Es ist 
05 (Sa):-0 = 0 (214) 
und 
8 ae (Stade (215) 
: a2 In den Abb. 5 und 6 ist S%, in Abhingigkeit von 
0 t* fiir verschiedene Werte von y aufgetragen (volle 
93 Linien.) Die Kurven wurden durch numerische Inte- 
gration aus (213) gewonnen}. 
Fiir den analog (206) ttber die Kreisflache der 
“ad Last gemittelten Wert (w,,); = 0 + der Anfangssetzung 
af (wm)no4 =x | 2rx[w(r)e-0 dr (216) 
a? 7 ee 
e ‘ 
0 QW Ql GB 00! 0S que Gur 008 G09 G7 erhalt man mit (212) aus (165) 
Abb. 5. Zeitlicher Verlauf der Nachsetzungen S;* 8a 
(volle Linien) und §* (strichlierte inten), (Wm)e—0 4+ = 3n (l+ )qqa. (217) 


_} Die zahlenmaBigen Auswertungen wurden vom Mathematischen Labor der Technischen Hochschule in 
Wien durchgefiihrt, wofiir wir seinem Vorstand, Herrn Professor Dr RK. Inzinger und dessen Mitarbeitern 
bestens danken, : 
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Benutzt man zum Dimensionslosmachen von (w):=04 denselben Faktor wie in (212), so erhalt 
man aus (217) 


ana 218 
lle Aang BENS Ti (218) 


Man erkennt aus (214), (215) und (218), da® der Wert der mittleren Anfangssetzung unter der 
Kreisflache stets gréRer ist als der Endwert der darauffolgenden mittleren Nachsetzung. 


Sia 
1 ae Soa 
A 
“a SZ zs = | 
=o 
cps = 
a 
al cE Z a = 7 
Za ae 
LZ LE 
O74, 
V=USAA a 
Z Y 
eZ 7) 
V=02 GH: 
06 Z C4 
ZRY Bhs 
LA 7 Pv 
05 Lip pps | VL ae 
ZF YH 
AZBAGL 
A OD a 
CA LN 
Zp hheas 
UA Ma = 4 Bins 
7 
Z 
Ze 
43 os 
“01 4 G3 = O04 05 06 0708097 2 3 fa O 9 10 20 30 40 50 6070 89000 


Abb. 6, Fortsetzung der Kurven von Abb. 5, Zeitachse logarithmisch geteilt. 


b) Setzung unter einem starren Stempel. Zum Vergleich soll noch die Setzung unter 
einem starren, durchlassigen Stempel vom Halbmesser a naherungsweise berechnet werden. In 
diesem Fall ist in der Randbedingung (89) q(r, t) unbekannt und so zu bestimmen, daB 


rst) =.0° far. a (219) 
und 
(w).=0 =f(t) (220) 
r<a 
mit der Nebenbedingung 
f2 Regt, ta Ps (221) 
0 


wobei f(t) eine noch zu suchende Zeitfunktion und P die vorgegebene Gesamtlast ist. P sei zeitlich 
konstant. 


Bei Beqehtine von (133) lautet (219) im Laplaceschen Unterraum 


ice) 
Q(r,s) =f QP, 8) pJo(rp) dp =0 fiir r >a. (222) 
0 
Die Laplace-Transformation von (220) liefert 
(W).=0 = F(s) (223) 
r<a 
mit 
ee) 
F(s) = fe—*f(t) de. (224) 
0 
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Aus (138) und (223) ergibt sich 
e 
Op, 8) Sale) dp ie fan eo eee 
i —2») p(p — Jp? + %s) + (lL—») xs pS ad —W)ns 


0 


Gleichung (221) lautet im Laplaceschen Unterraum 


29 | r QO(r,is) dr = = ; (226) 
0 


Driickt man Q(r, s) gema® (133) durch Q(p, s) aus, dann erhalt man nach Vertauschung der Inte- 
grationsreihenfolge bei Beachtung von (207) die Beziehung 


| (ps) Ji(ap) dp = =. (227) 


Die Gleichungen (222) und (225) stellen zwei duale Integralgleichungen fiir die unbekannte Funk- 
tion Q(p, s) dar, wobei F(s) durch die Nebenbedingung (227) bestimmt wird. Da diese Gleichungen 
in allgemeiner Form nur durch unibersichtliche Reihenentwicklungen léshar sind}, soll eine 
Naherungsliésung angegeben werden. Diese stitzt sich darauf, daB, wie zunachst gezeigt wird, die 
Lésungen dieser Gleichungen fiir t = 0+ und t = 00 auf dieselbe Druckverteilung gq(r) fiihren. 
Diese soll dann naherungsweise fiir alle ¢ als giiltig angesehen werden. 

Fiir die Zeiten t = 0+ bzw. t = oo gelten folgende Umkehrformeln fiir die Laplace-Transfor- 
mation ?: : 

Wenn lim g(t) existiert, so ist 


t>0 
lim g(t) = lim [s &{g(t)}], ~ (228) 
t>0 $00 
bzw. wenn lim g(t) existiert, so ist 
t+ oo 
lim g(t) = a [s L{g(t)}] . (229) 
t—>oo s—> 
Bezeichnet man demnach 
lim [s Q(p, s)] = qo(p) » (230) 
s—>0oO 
wobei gemaB (228) 
co 
go(P) = J Lalr)k=04 7 Jo(rp)ar (231) 


_ und 


lim [s F(s)] =f) 


dann erhalt man aus (225) nach Multiplikation mit s? und Ausfiihrung des Grenziiberganges s > 00 
die Beziehung 


| qo(P) Jo(r p) dp = aay firr <a. (233) 
0 


Nach Multiplikation mit s liefert derselbe Grenziibergang aus (222) 


Jap) pJlrp)dp=0 fir r >a. (234) 


Die Lésung des Systems (233) und (234) lautet 3 


aT ES ee (235) 


1 Vel. z. B. S. 111ff. des in FuBnote 1 von S. 236 zitierten Buches. 


2 G. Doetsch, Anleitung zum praktischen Gebrauch der Laplace-Transf i ii 
3 Vgl. S. 117 des in FuBnote 1 von Seite 236 zitierten Hachoal PSE Nari anrecaaceaine 
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Damit liefert die Umkehrung von (231) laut (93) 


20) | Oe | 
[g(r)]:=04. = Ser [ sin ap: J(r p) dp =) 2«(1+») ya? —72 Tea, (236) 
9 0 fur a. | 
_ Die Nebenbedingung (221) gibt 
P 
FO) > ree (1+). (237) 
Somit ist 
| a : fiir r< | 
1 
[a(r)-0y =) 20 Ry aa (238) 
|o fir r>a | 
Um [9(r)]:-. zu finden, bezeichnen wir 
lim [5 Q(P, 5)] = Go (P) (239) 
mit >: 
= co 
GooP) = J [4(r)}—co rJo(r p) dr (240) 
und 
oe [s F(s)] =f (co). (241) 
Dann erhalt man analog zu (233) und (234) aus (225) bzw. (222) durch die Grenziiberginge s —> 0 
co 
— f(~) ws 
i qoo(P) Jo(r p) dp = py ees) fiir r<a (242) 
6 
bzw. 
Oe. 
J Go(P) P Jor p)dp =0 fiir r >a. (243) 
Diese beiden Gleichungen haben dieselbe Form wie (234) und (235). Es ist demnach 
= f(e) sinap 
TAP) eae | (244) 
Die Riicktransformation gibt 
00 flo) 1 fa 
[q(r) }:—co == en | sin ap . Jer P) dp = | a (1—v?) a Yo —r? ur Tf =a a, | (245) 
. ‘2 0 fir r>a, | 
Die Nebenbedingung (221) fihrt auf 
1—).P 
f(co) = a (246) 
Semit ist 
P 1 re 
th. = (247) 


0 fir r >a. 
Es gilt also, wie gezeigt werden sollte, 
[a(r)]:=04+ = [9(7)]=0- (248) 
Aus (235) mit (237) bzw. aus (244) mit (246) folgt 


= a P sin 
Go(P) = GolP) =aaa (249) 


‘Wenn wir nun naherungsweise q(r) als zeitlich konstant ansehen, ergibt sich aus (249) 


Op. s) = 24G(p)} = = A. (250) 


2ua ps 
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Durch Einfiihren dieses Wertes in die Gleichung (138) bis (141) kann die zu dieser Druckverteilungen 
gehorige Losung ermittelt werden. Es soll hier jedoch nur noch die Setzung der Oberflache be- 
rechnet werden.! Die Gleichungen (138) und (250) ergeben fiir z = 0 


ee) 
~ P i 1 sin a p 
W,9=5— 401 =|?) | == Jy(r p) dp. 251 
o= ye HO ( ) | Gap = FES 5 Jolt p) dp (251) 
0 
Nach Rationalmachen des Nenners und Partialbruchzerlegung wird 
~~ ean cei } A 
W. iP BER (1 v) \P? ate OPS) p+ VP? +s Ju) dp : (252) 


z=0 =F xa—) | la - (1 — 27) p?+ (l—»)?x s 3 Dae 
0 


Mit Verwendung der fiir z = 0 spezialisierten Formeln (149), (151), (152) und (154) erhalt man 


co 
wih) 1— 2» pt Aes 
~ P pt v (i=) v p’t\|sin a p 
fe yy] [2 Sn eae ( ) erte (7/2) = J,(rp)dp. 
0 
(253) 
Daraus folgt die Anfangssetzung 
2 - 20 a (1 +9) fur ras 
Mah area At ala aes a 1 (254) 
f sage (1 +») are sin — fir r2a 
und der Endwert der Gesamtsetzung 
ii 
= a co 35% —*) fir r<a, 
= sina 
tw a (i =?) | AP. I(r p) dp = | (255) 
—.« (1 — v2) are sin — tit r22 as 
wa r 
Als ,,Nachsetzung“ soll wieder der Ausdruck 
S(r, t) = [w(r, t) — w(r, 0+)].—0 (256) 


bezeichnet werden. Mit (253) und (254) ist 


S(r, t) => «(1 | la—2y—0—9 Efe |/ 


1l— 2» p*t 


ie ie vy 1/p?t) | sin 
+ve Erfe ( =y - )| ar Jo(r p) dp. (257) 


Da fiir 0 <t < & die Setzung zufolge der gemachten Naherung auch fiir r <a noch von r ab- 


hangt, soll analog zu (206) wieder der Mittelwert S,,(t) tiber die Kreisflache eingefiihrt werden. 
Mit (207) wird 


$0 =4-0 +9 | a—29—a—» Bat = 


1—2»p*t pole 
5 (eae tare 2 rl 
tre OOF Bote (2 ENS ees dp. (258) 


1 Im folgenden werden die Grofen, die sich auf die Bel i i 
pyetiias sae ee 1e sich auf die Belastung durch einen starren Stempel beziehen, durch 
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Fuhrt man die Substitution a p = p, durch und verwendet wieder die dimensionslose Zeit nach 
(209), dann wird schlieBlich, wenn man statt p, wieder p schreibt 


— Sen. +9 | [a—2)—a— me (om) 


0 
iE 2 ze 
TiS reat Bete py) | Et) dp. (259) 
Daraus folgt 
co 
~ Pods. 1S i Pa(l La 
(Sp)ico = EE | REP I(p) dp = 72 OF OP) (260) 


0 


Nun lassen sich fiir die beiden behandelten Falle der konstanten Flachenlast unddes starren Stempels 
bei gleicher Gesamtlast 

Pi= an q (261) 
die mittleren Anfangssetzungen (w,,);.0, und (W,,):-0, bzw. die Endwerte der mittleren Nach- 
setzungen (S,,):;-< und (Sm):ioo. vergleichen. Aus (217) und (254) erhalt man 


(wm)t=0+ 32 

(@mi=04 322 eel (262) 
bzw. aus (211) und (260) 

(Sales 32 

aS oy Se || LS 263 

RG ae 372 ( ) 


Man erkennt daraus, da die beiden Grenzfalle des absolut nachgiebigen und des absolut starren 
Stempels bei gleicher Gesamtlast nahezu dieselben Werte fiir die tiber die Kreisflache gemittelten 
Anfangssetzungen und fiir die gemittelten Endwerte der Nachsetzungen liefern. Um auch den 
zeitlichen Verlauf vergleichen zu kénnen, beziehen wir analog wie friiher S,, auf (S,,):;-.. und 
erhalten fiir die dimensionslose Nachsetzung 


re Sm 4a Fa 
Sian Pa) — 20) Sm (264) 
Aus (259) ergibt sich 
Cc 
Ft Hleadalh Erfe ( /t*) 
~~" (1 — 24) P 
i 
| el see 
ey eae ear ye) [820 dp. (268) 


In Abb. 5 und 6 ist S* gemaB (265) in Abhangigkeit von der dimensionslosen Zeit t* nach (209) 
fiir verschiedene Werte von y dargestellt (strichlierte Linien). Man erkennt durch Vergleich mit 
den voll ausgezogenen Kurven fiir S;,, daB auch die zeitlichen Verlaufe der bezogenen mittleren 
Nachsetzungen fiir die beiden Grenzfalle des absolut nachgiebigen und des absolut starren Stempels 
nicht wesentlich voneinander abweichen. 


(Eingegangen am 19. Juni 1960.) 


Anschrift der Verfasser: 0. Prof. Dr. G. Heinrich und Doz. Dr. K. Desoyer, Wien 1V., Karlsplatz 13, 
Technische Hochschule. 
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Einige Auswerteverfahren bei spannungsoptischen Plattenuntersuchungen 
Von G. Haberland 


1. Einleitung. Mit Hilfe des Zweischichtenverfahrens oder des Reflexionsverfahrens lassen sich 
bei elastischen Platten spannungsoptisch die maximalen Torsionsmomente (d.h. die Differenzen 
der Hauptbiegungsmomente) und die Richtungen der Haupthiegungsmomente messen. Zur Be- 
stimmung der einzelnen Biegungsmomente kénnen rechnerische Auswerteverfahren herangezogen 
werden. Es wird gezeigt, daB sich praktisch alle fiir Scheiben bekannten Auswertemdglichkeiten 
auch bei Platten, die innerhalb der Kirchhoffschen Bedingungen auf Biegung beansprucht werden, 
ableiten lassen, sofern man — im Gegensatz zu den Scheiben — nicht von den Gleichgewichtshe- 
ziehungen, sondern von den Formanderungsbeziehungen ausgeht. In der Durchfiihrung der Aus- 
wertung besteht kein prinzipieller Unterschied zwischen Scheiben und Platten. 

Unter dem Begriff ,,Platte‘ verstehen wir einen ebenen Flachentrager, dessen Dicke im Vergleich 
zu den tibrigen Abmessungen des Kérpers klein ist und der einer Belastung unterworfen ist, die eine 
Kriimmung der urspriinglich ebenen Mittelflache hervorruft. 

Wir beschranken uns dabei auf solche Krafte, die in Richtung der Plattennormalen angreifen. 
Die Angriffspunkte dieser Krafte denkt man sich in die Plattenmittelflache verlegt. Denselben 
Korper bezeichnen wir dann als ,,Scheibe“‘, wenn er nur durch Kr§afte in seiner Ebene beansprucht 
wird und seine Mittelflache bei der Verformung eben bleibt. 


2. Grundlagen. Fiir ein rechtwinkliges Plattenelement, welches so in einem rechtwinkligen 
Koordinatensystem (x, y, z) liegen mége, daB die unverformte Plattenmittelflache und die (x, y)- 
Ebene zusammenfallen, gelten folgende drei Gleichgewichtsbeziehungen zwischen den Biegungs- 
momenten M,, M,, den Torsionsmomenten T’,,, den Querkraften Q,,, Q,, und der Belastungs- 
funktion! P(x, y): 

6 Txy 


aM, 
Ox = dy Uae (1) 
@My OTxy 
ay aac See (2) 
OQx x OQyz 
be. by oe (3) 


Weitere Gleichungen zur Bestimmung der fiinf unbekannten Funktionen M,,, M. y? Le Os Qys 
gewinnt man aus den Beziehungen, die zwischen den Spannungen und den Verformungen eines 
elastischen Kérpers bestehen (Hookesches Gesetz) sowie aus bestimmten vereinfachenden Annahmen 
tiber den Verformungszustand (Kirchhoffsche Bedingungen). Diese Annahmen bestehen darin, daB 
im Falle einer im Verhaltnis zu den Liangsabmessungen kleinen Plattenhéhe h und einer im Ver- 
gleich zu h kleinen Durchbiegung w diese als unabhangig von z angesehen wird und alle Punkte 
einer Normalen zur Plattenmittelflache auch nach der Verformung noch auf einer Geraden, der 
Normalen zur verformten Mittelebene, liegen. Das bedeutet eine Vernachliassigung der durch die 
Schubspannungen T,, und Ty, hervorgerufenen Verformungen. 

Fir die Verschiebungen u, v gelten somit die Beziehungen 


a ow ow 
u> z "Sco Py v=—Zz oy > (4) 
so das man schlieflich fiir die Verzerrungskomponenten folgende Ausdriicke erhalt : 
eS ow 
a : Pe ee ox? ’ 
eee ee 
y oy ay? ” (5) 
ee ov o?w 
Vey Sag V de arg 


1 Im allgemeinen ist es in der Plattentheorie iiblich, die Belastungsf i i i 
‘ en is : : gsfunktion mit p(x, y) zu bezeichnen. 
Wir verwenden dafiir die Bezeichnung P(x, y), da wir den Buchstaben Pp spater noch Sr parte verbrauchen, 


\ 
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Die Verzerrungskomponenten sind durch das Hookesche Gesetz mit den Normalspannungen o, 
und o, bzw. den Schubspannungen Ty verkniipft: 


Seni 1 2 (1 
Eg = = (Cx Hy) re aga Orv anll Ox)s Vey esta (6) 
Darin bedeuten yu die Poissonsche Konstante (u <0,5) und E den Elastizitatsmodul. Die Normal- 
spannung o, wurde neben o, und o, vernachlassigt. 
Fir die Spannungen folgen aus (5) und (6) die Beziehungen 


aes E . 02w 02w 

% 1 — ax! oy? 2 

sees E s 02 w C2 w 

y lp” Vo ! oe}? (7) 
ee. E ‘ 6?w 

ea lt+pw° 6x dy” 


Die auf die Lingeneinheit bezogenen Biegungsmomente M,, M y und die Torsionsmomente T,, 7 
ergeben sich nun auf folgende Weise: 


+h/2 
02w 02w 
M, = | o,2d = —B( ze +4 Ge) 
—h/2 
+h/2 
aw 02w 
M, = [ eee= B(T Hu oa) (8) 
—h/2 
+h/2 
62 w 
Ty = EL Sard Arf) 01s re 
—hj2 
Darin bedeutet 
E hs 
B= aa) 


- die Biegungssteifigkeit. 
Zwischen den 6rtlichen Ableitungen der Biegungsmomente und der Torsionsmomente besteht 
- danach bei Platten mit konstanter Biegungssteifigkeit folgender Zusammenhang: 


8(Mz —p My) _ OT xy 
2 EM) =a + u) | 0) 
und 
Fihren wir noch die Ersatzmomente m, und m, ein, die folgendermafien definiert werden: 
e2 
m, = —(1—y*) B—; =M,—pM,, (11) 
2 6?w es 12 
dann kénnen wir an Stelle von (9) und (10) schreiben 
=a Omx me OTxy 
ee te (13) 
omy oT xy ’ 


Diese Gleichungen werden allen unseren weiteren Betrachtungen, die die vollstandige Bestimmung 
des Biegespannungszustandes aus spannungsoptisch bestimmbaren GréBen zum Ziele haben, zu- 
_grundegelegt. In der Form entsprechen diese Gleichungen den Gleichgewichtsbeziehungen des 
-ebenen Spannungszustandes und werden auch entsprechend gehandhabt, ihrem Inhalt nach sind 
sie aber wesentlich von jenen verschieden. Um dies deutlich werden zu lassen, wurden hier noch 
_einmal die Grundlagen der klassischen Plattentheorie kurz aufgezeigt. 


EN 
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Wahrend die Gleichgewichtsbeziehungen des ebenen Spannungszustandes wie auch die hier 
angegebenen Plattengleichgewichtsbeziehungen (1) bis (3) — abgesehen von der Vernachlassigung 
kleiner GréBen héherer Ordnung — immer streng gelten, liegen den Gleichungen (13) und (14) die 
in der klassischen Plattentheorie iiblichen Annahmen itiber den Verformungszustand zugrunde. 
Die Méglichkeit, die ersten Ableitungen der Biegungs- und Torsionsmomente miteinander zu ver- 
kniipfen, beruht auf den Gleichungen (4); denn aus ihnen folgt, daB man fiir die Verzerrungen bei 
Platten schreiben kann 

Gee _ 1 Yay Oty 1 Opny (15) 
oy eee ax yeas 


Ohne die Gleichungen (4) ware erst eine Verkniipfung der zweiten Ableitungen der Verzerrungen 
méglich, was zu keiner fiir spannungsoptische Zwecke brauchbaren Beziehung fiihrt. 

Ebenso kénnen wir die Gleichgewichtsbeziehungen (1) und (2) nicht unmittelbar benutzen, da 
sie ja die unbekannten Querkrafte enthalten, die sich durch Einsetzen der Gleichungen (8) erst 
formelmaBig aus diesen Gleichgewichtsbezichungen ergeben. 

Wir diirfen uns also nicht dariiber hinwegtauschen, dafB bei den hier beschriebenen spannungs- 
optischen Untersuchungsméglichkeiten elastischer Platten die gleichen — streng genommen unzu- 
langlichen — Voraussetzungen enthalten sind, die auch einer theoretischen Behandlung nach der 
klassischen Kirchhoffschen Theorie zugrundeliegen. Das ist aber auch bei vielen anderen experi- 
mentellen Verfahren der Fall, z. B. dann, wenn die Biegungsmomente aus Kriimmungsmessungen 
nach (8) ermittelt werden. 

Der Zusammenhang zwischen den Momenten verschiedener Schnittrichtungen wird in bekannter 
Weise durch den Mohrschen Momentenkreis veranschaulicht. Insbesondere stehen die Biegungs- 
momente M,, M, sowie das Torsionsmoment T',, in folgendem Zusammenhang mit den Haupt- 
biegungsmomenten M,, M,: 


M, = M, cos* » + M, sin* ¢ , 


ci Pe M, sin? pg + M, cos? , (16) 
(M,—M,). 
oe sin2@, 


worin mit » der Winkel gemeint ist, den die Richtung des Hauptbiegungsmomentes M, mit der 
x-Achse bildet. 


1 : . 
ra (M, — M,) bedeutet das maximale Torsionsmoment T,,,,,. Fiihren wir auBerdem 


M,-+ M,=p und» M,— My =2 Ta = 4 (17) 
ein, dann lassen sich die Gleichungen (16) auch folgenderma8en schreiben: 
Biel Oped 
M,=> +> cos29, | 
19 q 
My= = = 2082. (18) 


T,y= Tmax Sin 2 ~ = 4 sin 2 P| 


Fiir die Ersatzmomente m,, m, in (11) und (12) kénnen wir mit (18) schreiben 


m, = (1—p) + + (1+) 4 cos2 9, 


(19) 
m, =(1—p) + —(1 +p) {cos 2. 
Setzen wir noch (19) und die dritte Gleichung (18) in (13) und (14) ein, dann ergibt sich 
CPs Wee to quar 0q op op. 
Tel = E sin 2 p — a cos2y+2q te cos2y + a sin2¢)], (20) 


Opae ase oq OG ox 6 r ta] 
Ge — Tap [an O82 P+ Gy sin 2 + 2g(— FE cindy + con?) (21) 
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3. Die spannungsoptischen Mefgréfen. Da dic Biegungsspannungen nach dem Geradliniengesetz 
[vgl. (7)] ttber die Plattenhéhe verteilt sind, kann eine homogene Platte bei einer Beobachtung im 
polarisierten Licht keinen spannungsoptischen Effekt zeigen. Die Spannungen haben zu beiden 
Seiten der neutralen Flache verschiedene Vorzeichen, so daB die optischen Gangunterschiede in der 
Summe Null ergeben. Man erhalt aber ein spannungsoptisch beobachtbares Bild, wenn man die 
Modellplatte aus zwei Schichten verschiedenen Materials mit unterschiedlicher spannungsoptischer 
Wirksamkeit zusammenklebt (Zweischichtenverfahren) oder aber in die Plattenmittelflache eine 
reflektierende Schicht einklebt und die Platte im reflektierten Licht beobachtet (Reflexionsver- 
fahren)!. In beiden Fallen erhalt man Isochromaten, aus denen man die Tmax-Werte ermittelt, und 
Isoklinen, die die Richtungen der Haupthiegungsmomente angeben. Nach derletzten Gleichung (18) 
werden also die Torsionsmomente T,, iiber die ganze Platte experimentell bestimmt. 

Die Theorie der Zweischichtenplatte ist bereits in der Literatur ausfiihrlich beschrieben worden, 
so da8 hier nicht naher darauf eingegangen zu werden braucht®*%56, Es soll hier der Hinweis 
gentigen, da wir fiir die Zweischichtenplatte lediglich resultierende elastische Konstanten anzu- 
nehmen haben, die sich aus den elastischen Konstanten und Dicken der einzelnen Schichten her- 
leiten lassen, falls diese unterschiedlich sind. Die ermittelten Biegungs- und Torsionsmomente 
diirfen wir als einer homogenen Platte zugehérig betrachten (wodurch das Zweischichtenverfahren 
letzten Endes erst sinnvoll wird). Favre, Gilg und Schumann haben iiber verschiedene spannungs- 
optische Untersuchungen mit Hilfe des Zweischichtenverfahrens berichtet” *®®, Da von ihnen die 
einzelnen Haupthiegungsmomente direkt mit Hilfe eines Interferometers gemessen wurden, er- 
ubrigte sich bei ihnen die Anwendung eines Auswerteverfahrens zur Bestimmung der einzelnen 
Biegungsmomente. Kuske und Hirschfeld haben vorgeschlagen, Zweischichtenplatten als Ganzes 
zu beobachten, d.h. Isoklinen und Isochromaten aufzunehmen, so wie es bei Scheibenversuchen 
allgemein itiblich ist! 1, und was in der experimentellen Durchfihrung auSerdem einfacher ist. Da 
aber hierbei nur die maximalen Torsionsmomente und die Richtungen der Hauptbiegungsmomente 
gemessen werden kénnen, sind rechnerische Verfahren notwendig, um aus den MeSgréBen die ein- 
zelnen Biegungsmomente iiber die ganze Platte zu bestimmen. Der Wiedergabe einer Reihe der- 
artiger Auswerteverfahren sind die folgenden Abschnitte gewidmet. Die erforderlichen Grundlagen 
dazu sind-bereits im vorigen Abschnitt entwickelt worden. 


4. Die Bestimmung der einzelnen Biegungsmomente lings geradliniger Wege. Diese Auswertung 
beruht auf einer unmittelbaren Anwendung der Gleichungen (13) und (14). Sie entspricht in ihrer 
Durchfihrung dem bei Scheiben iiblichen Schubspannungsdifferenzenverfahren, welches auf den 
Gleichgewichtsbeziehungen des ebenen Spannungszustandes aufbaut. Sollen die Biegungsmomente 
itiber die ganze Platte bestimmt werden, dann ist es zweckmafig, die Torsionsmomente T,, fiir 
hinreichend viele Schnitte x = konst. bzw. y = konst. aufzutragen und diese Kurven graphisch 
zu differentiieren. Die Bestimmung der Ersatzmomente m, und m, im Innern der Platte erfolgt 
dann durch Integration der értlichen Anderungen der Torsionsmomente lings geradliniger Wege 
von Randpunkten P, bzw. P; aus, in denen die Ersatzmomente aus den jeweiligen Randbedingungen 
bekannt sind: 


Te 

(m)e = (L+H) | dy + (mdpo (22) 
P, 
fz) 

(mye = (L+H) | GP de + (m)es- (23) 
Py 


1 J. N. Goodier u. G. H. Lee, J. Appl. Mech. 8 (1941) S. A 27. Dazu Diskussionsbeitrag von R. D. Mindlin, 
. Appl. Mech. 8 (1941) S. A 187. 
: Pg. Favre - B. Gilg, Schweizerische Bauztg. 68 (1950) S. 253 und S. 265. : ¢~< 
8 B. Gilg, Experimentelle und theoretische Untersuchungen an diinnen Platten. Diss. E. T. H. Ziirich 
1952 (Publication No. 5 du Laboratoire de Photoélasticité de I’Kcole Polytechnique Fédérale, Ziirich 1952). 
B. UW. Kopoues, UnKxenepubiii COopHuk 22 (1955), S. 98. 
G. Haberland, Przeglad Mechaniczny 19 (1960) S. 120 und S. 153. 
H. Schwieger u. G. Haberland, Z. angew. Math. Mech. 39 (1959) S. 85. 
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: e Dos ae Gihiteanti Etude expérimentale de la répartition des moments dans une plaque oblique 
fléchie en fonction de l’angle par les cétés. Bulletin Nr. 9 de la Société Frangaise des Mécaniciens (1953). 
10 A, Kuske, Z. VDI 94 (1952) S. 745. 
11 K. Hirschfeld, Bauingenieur 25 (1950) S. 455. 
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Bei der praktischen Anwendung dieser Gleichungen kann es mitunter auch angebracht sein, die 
Reihenfolge von Differentiation und Integration zu vertauschen; denn wir kénnen ja 0/0x baw. 0/dy 
auch vor das Integral setzen. Es kann also auch zuerst tiber die Torsionsmomente T,y in der einen 
Richtung integriert werden, um erst hernach diese Integralwerte in der zur Integrationsrichtung 
senkrechten Richtung zu differentiieren. 

Soll eine Platte vollstandig ausgewertet werden, dann ist es zweckmafig, zunachst einen Schich- 
tenplan fiir die Torsionsmomente T,,, zu zeichnen, was auch graphisch erfolgen kann mit Hilfe des 
Mohrschen Kreises entsprechend der Konstruktion des ,,Schubspannungshiigels‘‘ bei Scheibenver- 
suchen, wie es Kuske vorgeschlagen hat’. 

Es geniigt, ein Ersatzmoment, z. B. my, durch Integration zu bestimmen. Man erhalt dann m, 
gem4B (19) und (17) aus der Beziehung 


m, = my +2 (1 +4) Tax cos 2 yp. (24) 
Die Biegungsmomente M, und M, erhalt man folgendermafien aus den Ersatzmomenten [vgl. (11) 
und (12)]: 
1 
IS pope TH M,), 
5 (25) 


Dieses Verfahren wurde zuerst von Schwieger vorgeschlagen®. Es hat sich inzwischen bei ver- 
schiedenen Plattenuntersuchungen bewdhrt**. Es wird wohl auch dasjenige Verfahren sein, dem 
vor den anderen Auswerteméglichkeiten in vielen Fallen der Vorzug zu geben sein wird. 

Dieses Verfahren laBt sich auch leicht auf Platten mit veranderlicher Dicke erweitern. Die Ver- 
kniipfung der Ersatzmomente m, und m, mit dem Torsionsmoment T,, fiihrt hierbei, wie man aus 
(8), (11) und (12) leicht findet, zu den Ausdriicken: 


ets ee 
m(B)=0 +4 a (a) 


Die Ersatzmomente im Punkte P ergeben sich wieder durch Integration von den Randpunkten P, 
bzw. P) aus: 


(26) 


(m,)p = (1 +4) Bp fe (F) & + Be (F), 
= (27) 
(m,)p =(1 +p) Bp | at a dx + Bp ls . 


Pp! 
0 
Als Anwendung wurden die Biegungsmomente einer Kragplatte mit linearer veranderlicher Dicke 


bestimmt®. 


5. Die Bestimmung der Momentensumme lings heliebiger Integrationswege. Zur Bestimmung 
der Momentensumme durch Integration langs eines beliebigen Weges sind im Prinzip schon die 
Gleichungen (20) und (21) geeignet. Da die Orientierung der (x, y)-Achsen in der Plattenmittelflache 
zunachst beliebig gewahlt werden kann, miissen gleichartige Beziehungen entstehen, wenn wir an 
Stelle von x, y das orthogonale, aber sonst beliebig in der Plattenmittelflache liegende System é, n 
einfiihren, wobei dann an Stelle von @ der beliebig veranderliche Winkel 6 tritt, den die Richtun gé 
mit der Richtung des Hauptbiegungsmomentes M, bildet. Wir wollen die Gleichungen (20) und (21) 
in diesem Sinne umschreiben, da es fiir die Anschauung bequemer ist, das einmal festgelegte Koor- 
dinatensystem x, y feststehen zu lassen. 


: A Kuske, Forsch. auf dem Gebiete des Ingenieurwesens 18 (1952) S. 113. 
E Bee - Schwieger, Bauplanung und Bautechnik 8 (1954) S. 174 und Wiss. Z. Univ. Halle, Math. Nat. 4 (1954) 
aes H. Schwieger u. G. Haberland, Bauplani und Bautechnik 8 (1 i i 
made art ioe Saas planung utec. (1954) S. 358 und Wiss. Z. Univ. Halle, 
4 Siehe FuBnote 5 von Seite 257. 
5 Siehe FuBnote 6 von Seite 257. 
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Ks seien s, und sy die Richtungen der Hauptbiegun ie mi 
: gsmomente M, und M,, die mit dem x, y-System 
den Winkel P bilden, # der beliebig verinderliche Winkel, den ide Ricken & mit Ac etches 
bildet, 7 sei die zu  senkrechte Richtung. Der Winkel 6 zwischen £ und Ss, ist dann (Abb. 1) 


6=0— 9. (28) 


Abb. 1. Bezeichnungen zur Integration lings eines beliebigen Weges. 


Wir beziehen zunachst die Gleichungen (13) und (14) auf die Richtungen!? é, 7. Es ist 


Omx dmx OF | Omx On dmx .. mx 
ee ay om ery: sin } + en cos? , ] 
omy a omy 6& omy on at § omy omy ° 
em ai, bx, Eo 
29) 
OTxy _ OT xy oé OTxy an _ OT xy OTxy . ( 
Reeeieates.| > 6 62 ee Sey 
OT xy , OTxy 0& OTxy Hy) a OT xy . OTxy 
terest packs exes ee sin O + bn cos? . 
Setzen wir diese Ausdriicke in (13) und (14) ein, dann ergibt sich 
Omz . omx oe 6Txy clay sas 
5 sin } + 7 cos #— (1 + yp) ( 5 cos 0 én sin | ==105 (30) 
omy ___ omy ‘ gehts OTxy « OTxy 250 
pe 0088 5 sin ® — (1 +4) ( RE sin? + mn cos 8) = 0. (31) 


Setzen wir fiir die Ersatzmomente m,, m, die Gleichungen (19) und fiir die Torsionsmomente T, , die 
dritte Gleichung (18) ein, dann erhalten wir zwei Beziehungen, von denen wir einmal die erste mit 
sin @, die zweite mit cos ? multiplizieren und addieren. Zum andern multiplizieren wir die erste 
Gleichung mit cos # und die zweite mit — sin } und addieren ebenfalls. Unter Beriicksichtigung 
von (28) erhalten wir als Ergebnis die beiden folgenden Gleichungen: 


op 1+ pu | éq Comes op. op : : 

ey E cond — HE sin 2d + 2a (Zp sin25 + 5 00828)], (32) 

Op lb pels ged ecg op _ & .. 
TRE e Zz sin 2 0 on cos26 +24 (Ze cos 20 on sin 24)]. (33) 


Auf den rechten Seiten dieser Ausdriicke stehen nur spannungsoptisch bestimmbare GréBen. Diese 
Gleichungen sind gleichbedeutend mit einem Hineindrehen der Richtungen x, y in die beliebigen 
Richtungen é, 7; denn sie gehen aus (20) und (21) mit dx > dé, dy > dy und y = — 6 hervor. 

Ist dg/0E = 0, so fallt die Richtung ¢ mit einer Isochromate zusammen. Im Falle dg/of = 0 
fallt die Richtung & mit einer Isokline zusammen. 


1 Hier sowie in den folgenden Abschnitten beschranken wir uns auf Platten mit konstanter Dicke und 
gehen deshalb von den Gleichungen (13) und (14) aus. Handelt es sich um Platten mit veranderlicher Dicke, 
d. h. mit ortsveranderlicher Biegungssteifigkeit B, dann hat man bei der Ableitung der einzelnen--Auswerte- 
methoden von (26) auszugehen. 
f 18* 
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Mitunter mag es zweckmaBig sein, die Ableitungen von q und @ nach s, bzw. s, einzufuhren, was 
wir noch tun wollen. Es ist 
(a) (a) a) : 
t= cos 6 + = sind, ] 
2 


0& 0s, 
a = “ sin 0 + ne coso , 
U] Sy 2 (34) 
ep ey op 
BE ba rapes re sind, 
op op oP 
AE ee — cos 6 
on ae 
Setzen wir diese Ausdriicke in (32) ein, dann erhalten wir 
ett | Pe Ot aj 2 See ) YF eceas 35 
ee” dl — yp | és, i 89 se ale as, ents 85 is) 


Wir integrieren diese Gleichung zwischen den Punkten 0 und | und erhalten 
1 


1 
| bs coso a sin ) dé + 2 | q ce sin 6 + =. cos 4) 7 (36) 
O85 5 Os, O85 
0 


Os, 


l+y 
1—wu 


Piz Pos 


oder 


dé. (36a) 


Pie Po = - | [sin 6 ( : + 2q wn) +0088 (5 +24 30 | 
4 | 


Das analoge Auswerteverfahren fiir Scheiben wurde von Foéppl angegeben!. 

Auch die Einfiihrung der Ableitungen von q und nach irgendwelchen anderen ausgezeichneten 
Richtungen in (32) und (33) kann mitunter angebracht sein [vgl. dazu die spateren Gleichungen (95) 
und (96)]. 


6. Die Integration lings einer Haupthiegungsmomentenlinie. Bei Scheiben ist auch das Ver- 
fahren gebrauchlich, die einzelnen Hauptspannungen durch Integration langs der Hauptspannungs- 
linien zu ermitteln. Dazu geht man von den auf die Hauptspannungsrichtungen bezogenen Gleich- 
gewichtsbeziehungen aus. Der diesem Verfahren analogen Auswertung bei Platten wollen wir uns 
jetzt zuwenden. Die dazuerforderlichen Beziehungen gehen als Spezialfalle aus (32) und (33) hervor. 

Es falle € in die s,-Richtung, 7 in die s,-Richtung. Setzen wir also 


df» ds,, dy— ds, .0=0 
in (32) und (33) ein, dann ergibt sich 


op itp /a ae \ 
és, 1—yp (et 2456) (37) 
ép _1l+p/__ oq op 
Soe ga +245) (28) 
Setzen wir gemaB (17) fiir p wieder die Momentensumme, fiir q ihre Differenz ein, dann wird aus (37) 
P a 
a; (Ma —» My) = (1 + 4) (M, — M,) 52 89) 
und aus (38) 
a @ | 
ge, (Mi —# M,) =(1 +») (My — Mh) 5. (40) a 
Mit Einfiihrung der Hauptersatzmomente 
: m, = M,—y M,, m, = M,—yp M, (41) 
wird daraus 
om 6 
$5, = +H) (My — My) 3 (42) 
und 
om, 7) 
35, = (L+H) (Ma — My) 5 - (43) 


1 L. Foppl, Z. angew. Math. Mech. 34 (1954) S, 454. 


am 
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Fiir die praktische Verwendung ist es mitunter vorteilhaft, diese Gleichungen in Differenzenform 
zu schreiben: 


Am, = (1 +) (M, — M,) 222 , (44) 
As, 
Am, = (1 +4) (M,— M,) 222 (45) 


Die Bedeutung der einzelnen GréBen ist aus Abb. 2 ersichtlich. Gleichung (44) gilt fiir ein Fort- 
schreiten von A nach B, Gleichung (45) fiir ein Fortschreiten von A nach C. Als Ausgang der 
Integration hat man einen Plattenpunkt zu wahlen, in dem M, — M, bekannt ist. Dieser Punkt 
wird gewohnlich am Plattenrande liegen. 


Lsokline p 
Abb. 2, Dreieckiges Plattenelement, gebildet aus Isokline Abb. 3. Trajektorien der Hauptbiegungsmomente bei 
und Hauptbiegungsmomentenlinien, einer rotationssymmetrisch belasteten Kreisplatte, 


7. Die Auswertung bei rotationssymmetrisch belasteten und gelagerten Kreisplatten. Im Falle 
einer zentrisch oder rotationssymmetrisch belasteten und lings des Umfanges gleichférmig ge- 
stiitzten Kreisplatte mit dem Halbmesser a sind die Haupthiegungsmomentenlinien die konzen- 
trischen Kreise um den Plattenmittelpunkt und die Radialstrahlen. Wir bezeichnen die Richtungen 
der Haupthiegungsmomente M, = M, mit s,, die der Haupthiegungsmomente M, = M, mit s, 
(Abb. 3). Fiir diesen Fall geht die Beziehung (42) wegen 

foie te es ealeoagh dé "Ar (46) 


Os5 On r 
(0. ist der Kriimmungshalbmesser der s,-Linie) iiber in 


dm, — M),) 


us (My 
ic Pe (pa (47) 
Durch Integration vom Rande aus folgt 
| M,—M 
(M,—1 Myo— 1 = (ms) = (1 +1) | Mar + (my), (48) 


Eine Differentiation gemessener GréBen entfallt somit in diesem Falle, da die Kriimmungen der 
Hauptbiegungsmomentenlinien durch die Rotationssymmetrie gegeben sind. _ 
Die einzelnen Haupthiegungsmomente ergeben sich nun leicht folgendermafen: 


M,=M, =72=4 und M, =M, =" mit q=M,—M,. (49) 


Im Plattenmittelpunkt ist M, = M,, d.h. q = 0. Aus (48) und (49) ergibt sich fiir die dortigen 


Haupthiegungsmomente 
0 


(m,), _ 9 1t+typ| M,—M, (m2), 4 
(M,), =0 = (M,),; =0 = yr plea Sem Ei 9 (50) 


wobei zu beachten ist, daB (M, — M,)/r selbst in diesem Punkte unbestimmt wird. ; 
Der Randwert (m,),—. ergibt sich aus den Randbedingungen. Sie lauten fiir einen frei drehbar 


aufliezgenden Rand 
e p10 und (M)). =, ==0 (51) 
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und fiir einen eingespannten Rand 


W,..=0 und ae =0. (52) 
Fiir die Biegungsmomente einer rotationssymmetrisch belasteten Kreisplatte gelten die Ausdriicke 
M, = M,=—B(@ te M, = M, = Blues cee) (53) 
Das Ersatzmoment m, wird also 
m, = M,—w M,=—BQ—,)—™. (54) 
Fiir den frei aufliegenden Rand erhalten wir wegen (M,),_ 4 = 0 
(75), — a = (Mg) =o = — Grae = — 2 (Tax) aa: (55) 
Am eingespannten Rand wird mit der Randbedingung (52) 
(ms), =a =9. (56) 


Auf die richtige Wahl der Vorzeichen ist besonders zu achten, da die Differenz der Haupt- 
biegungsmomente nur absolut gemessen wird. Von dieser Auswertemdglichkeit haben Schwieger 
und Reimann bei der spannungsoptischen Untersuchung einer quergestoBenen Kreisplatte Gebrauch 
gemacht'. 


8. Die Integration langs einer Haupttorsionsmomentenlinie. Ein weiteres bei Scheiben tibliches 
Auswerteverfahren beruht auf der Bestimmung der Spannungssumme durch Integration langs der 
Hauptschubspannungslinien. Ganz analog laBt sich bei Platten die Momentensumme durch Inte- 
gration langs der Haupttorsionsmomentenlinien bestimmen, wie jetzt gezeigt werden soll. Bekannt- 
lich bilden die Richtungen der Haupttorsionsmomente mit den Richtungen der Hauptbiegungs- 
momente in jedem Plattenpunkte einen Winkel von 45°, so da8 man das orthogonale Trajektorien- 
netz der Haupttorsionsmomente sofort aus dem der Hauptbiegungsmomente gewinnen kann. 

Bezeichnen wir die Richtungen der Haupttorsionsmomente mit t, und t, und die zugehérigen 
Linienelemente mit dt, und dt, (Abb. 4). Wir gehen wieder von den allgemeinen Gleichungen (32) 
und (33) aus und lassen ¢ mit der t,-Richtung, 7 mit der t,-Richtung zusammenfallen, setzen also 


dé — dt, , dyn — dty , 6 = 45° 
in (32) und (33) ein. Es folgt 


ea pte 2 429%), (57) 

gp = pt e(_ 2 _ 9g M8) (58) 

Fiir eine praktische Verwendung schreibt man diese Gleichungen in folgender Differenzenform: 
Ap =1*#(— dq 48429 de), (59) 

dp = 74 (— 4,4 42-24 do). (60) 


Darin gilt A, fiir ein Fortschreiten von A nach B und A, fiir ein Fortschreiten von A nach C (Abb.5). 


Tsokline p 


Abb. 4, Die Richtungen t, und t, der Haupttorsionsmomente. Abb. 5. Dreieckiges Element, gebildet aus Isokline 


und Haupttorsionsmomentenlinien. 


1 H. Schwieger u. V. Reimann, Z. angew. Math. Mech. 39 (1959) S. 198. 
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Die Gleichungen (59) und (60) sind dann den Gleichungen (44) und (45) vorzuziehen, wenn die 
Isokline mit der s,-Richtung einen sehr kleinen Winkel bildet, so daB das Dreieck in Abb. 2 duBerst 
langgestreckt wird und umgekehrt. 

op ep op op 


Zwischen den GréBen ; 
oben Gt,’ bt,” Os, ds, besteht noch folgender Zusammenhang: 


OP ee CPR 45° a OP in 45° — : y2 (2 aca 
$) 


Ot, Os, Os Os 
op op op 1 é é (61) 
ees gs rs eae 28 eee 2 fet OP a ee 
Bt, ee eee 3 |e ( be, cr) 
bzw. 
eye 2). 
Os, 2 Cty. at, 
F) 1 F) 8 (02) 
Pa, > (oP , op 
és, 2 y2 a . zr) ‘ 


9. Ein weiteres Auswerteverfahren zur Bestimmung der Momentensumme. Das Verfahren zur 
Bestimmung der Momentensumme, das jetzt beschrieben werden soll, stellt das Analogon dar zu 
einem von L. Féppl fiir Scheiben angegebenen Verfahren. Fiir Scheiben sind dieses Verfahren 
sowie die Verfahren, die auf einer Integration in Rich- 
tung der Hauptspannungslinien bzw. der Hauptschub- 
spannungslinien beruhen, in zusammenfassenden Dar- 
stellungen iiber Spannungsoptik beschrieben! 2 % 4, 


Man geht davon aus, daB man den Winkel y, den die 
Isokline mit der Hauptbiegungsmomentenrichtung s, 
bildet, beiderseits von s, auftragt. Diese beiden neuen 
Richtungen werden mit s und ¢ bezeichnet, die ent- 
sprechenden Linienelemente mit ds und dt (Abb. 6). 


Wir multiplizieren nun (37) mit sin y> (38) mit cos y Abb. 6, Die Richtungen s und ¢ zur Bestimmung 
5 : der Momentensumme. 
und addieren beide: 


op . op site eq. eg l+uyu op : ap 
oe. tT os, cos y artes Bes geet y) sa = 2 4(% siny Brrr y). (63) 


Die linke Seite dieser Gleichung bedeutet aber 0p/ds, der erste Klammerausdruck der rechten Seite 

ist gleichbedeutend mit — 0q/0t. Der zweite Klammerausdruck schlieBlich mu verschwinden, da 

er die Anderung des Winkels @ langs der Isokline angibt, p aber langs der Isokline konstant ist. 
Wir erhalten somit 


ee On (64) 
as 1—p @t © 


Mit Hilfe der (s, t)-Linien, die unmittelbar aus dem Isoklinenfeld gewonnen werden (die t-Linien 
stehen iiberall senkrecht auf den Isoklinen), kann somit die Momentensumme in s-Richtung erhalten 
werden, wenn man die Anderung der Haupthiegungsmomentendifferenz in der t-Richtung ermittelt 
hat. 

Einige Sonderfalle sollen noch erwahnt werden. Falls die Isokline die Hauptrichtung s, in 
einem Punkte beriihrt oder mit ihr langs einer Linie zusammenfallt, also y = 0 wird (z. B. Symme- 
trieschnitt), fallen die Richtungen s und t mit der s,-Richtung zusammen. Es wird dann aus (64) 

C7) 


o2(M22 os ee 65 
Fe, (Ma — t Ma) =, = 9» (65) 


was wir auch aus (43) mit 0~/ds, = 0 ablesen kénnen. 


1G, Mesmer, Spannungsoptik. Berlin 1939. ‘ ; 

2 L. Féppl u. E. Mének, Praktische Spannungsoptik. Berlin-Géttingen-Heidelberg, 1. Aufl. 1950, 2. Aufl. 
1959. 

3 Siehe FuGBnote 1 von Seite 258. : ie 

4 Fine ausfihrliche Beschreibung der Integration langs einer Hauptspannungslinie bzw. Hauptschub- 
spannungslinie findet man auch in M. M. F. rocht, Photoelasticity, New York, Bd. iT: 1941, Bd. II: 1948 und 
in C. B. Biezeno u. R.Grammel, Technische Dynamik, Bd. 1, 2. Aufl. Berlin-Géttingen-Heidelberg 1953. 
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Fallt die Hauptrichtung s, mit der Isokline zusammen, wird also y = 90°, dann fallt die s-Rich- 
tung in die positive s,-Richtung, die t-Richtung in die negative s,-Richtung, und es folgt aus (64) 


og = ee i es 66 
Emer pM) eS. =O ou 
was wieder aus (42) mit 0g/ds, = 0 folgt. 
Aus der Gleichung der Isokline: 
oy .. op #¢ 
nee Parse oe (67) 


entnimmt man, daf fiir beliebige Werte von y die Beziehung 


op (OP _ Qe 
rears a / és. 1 fe: 
besteht, wobei 0, und 9, die Kriimmungshalbmesser der Haupthiegungsmomentenlinien bedeuten. 
Insbesondere fiir y = 45° wird 0, =@. Ist, = Cound y + 0°, dann ist auch Q, = Cound umge- 
kehrt fiir 9, = co und y + 90° ist auch 0, = 00. Ist der Kriimmungshalbmesser einer Hauptbie- 
gungsmomentenlinie unendlich, dann verlauft sie entweder geradlinig oder sie besitzt einen Wende- 
punkt. 


10. Die Anwendung des Neuberschen Auswerteverfahrens bei Platten. Neuber hat ein Verfahren 
zur Konstruktion der Linien konstanter Spannungssumme bei Scheiben beschrieben!, das sich auch 
fiir die Konstruktion der Linien konstanter Momentensumme bei Platten erweitern la8t. Man geht 
dazu von vier orthogonalen Netzen aus (Abb. 7). Es bedeuten s,, s, wieder die Haupthiegungsmo- 
mentenlinien, s, sind die Linien konstanter Momentensumme p, s, ist die dazu senkrechte Kurven- 
schar, s; sind die Linien q = konst. (Isochromaten), s, wieder die dazu senkrechten Linien, und 
schlieBlich s, und s, sind die Isoklinen und die zu ihnen 
orthogonalen Kurven. Die Bezeichnung der Winkel 
zwischen s, und den iibrigen Richtungen ist aus Abb. 7 


ersichtlich. 
Nun wird 
op op OD as 6 : 
PP cosa — SP sina =— sinx (69) 
Os, O83 Os, OS, 
und 
ép Op me op (7) 
= sina + 2 cosa =~ cosa (70) 
G85 OS. OS, CS, 
Abb. 7. Die Bezeichnungen fiir die Tangentenrichtungen 
und -winkel bei den vier orthogonalen Netzen zur An-' Wepre 
wendung des Neuberschen Verfahrens. ép 


Os, 
Der Ausdruck dp/s, stellt den als positiv angenommenen Gradienten von p dar. Beim Ubergang 


von einer Linie p = konst. zur nachsten nehme p immer um den konstanten Betrag m, zu. Ist 
a der Abstand zwischen den Linien p = konst. und p + m, = konst., so gilt naherungsweise 


po Se 
; OS, Fe aa (71) 
Es ist somit 
op Oo: op ™o 
oo —— se —- = — : 
36, ina , a a (72) 
Nun ist definitionsgema8 auch 
oq op 
Oss —— und Os; > 0 > 
so da8 wir entsprechend (72) auch schreiben kénnen 
og ™ |. aq __ ™ 
PIN aes font foie co ce (73) 
und 
oe = — ma sin sa a 
Pa a Vo 08, = ia cos Y > (74) 


Lee Féppl u. H. Neuber, Festigkeitslehre mittels Spannungsoptik. Miinchen und Berlin 1935. 
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wobei q beim Ubergang von einer Linie g = konst. zur nachsten wieder um den konstanten Be- 
trag my) und y beim Ubergang von einer Isokline zur nachsten um den Betrag ¢ zunehmen sollen. 
Der Abstand zweier q-Linien ist mit b, der Abstand zweier Isoklinen mit c bezeichnet. Es werden 
also Hohenschichtlinien von p und q betrachtet, auf denen p und q ganzzablige vielfache Werte 
von my, annehmen, und Isoklinen, langs denen ¢ ein ganzzahliges Vielfaches von ¢ betragt. 


Die Ausdriicke (72), (73), (74) setzen wir in (37) und (38) ein und erhalten 


™M . Sse ™M . E 
— sin hallf Fe sin B +24 cosy), 
(75) 
1 
** cosa = 1TH / cos 6 —2q : sin 7] 
oder nach Multiplikation mit ( PE ] ae 
1+ pu] m 
te ibe é eb 
—(=4) ; ee ag ae OSE ly ag coe, 
: (76) 
ee be Ue 
ae Z cos % = —ccos B — 2 15. ey: 
Mit den Abkirzungen 
eSNG anew eb , 
(TH) be und ose (77) 
wird daraus 
a’ sina = csin§ —b’ cosy, a’ cosa~ = —ccosh—b'siny. (78) 


Das sind zwei Gleichungen zur Bestimmung der beiden Unbekannten a’ und «. Werden beide 
Gleichungen durcheinander dividiert, so erhalt man 


b' cosy —c sinB 


1B nt Brahe ENE BS (79) 
- Durch Quadrieren und Addieren der beiden Gleichungen (78) folgt 
a’? = b? + c? —2 b’ ccos(y + 90° — A). (80) 


Diese beiden Gleichungen haben die gleiche Form wie die von Neuber fiir Scheiben abgeleiteten 
Beziehungen. Die Bestimmung der p-Linien nach (79) und (80) erfolgt auf graphischem Wege, 
und zwar ganz analog wie bei Scheiben, so dafs auf das Buch von Féppl und Neuber verwiesen 
werden kann.1 

Bei Scheiben kommt der zu den Hohenschichtlinien der Spannungssumme orthogonalen Kurven- 
schar eine physikalische Bedeutung zu. Langs einer derartigen Linie ist nimlich dort die soge- 


nannte Drehung 
I /év Ou 
Grete be 


konstant? (q, ist die z-Komponente eines axialen Vektors @, der eine Drehung hervorruft, welche 
den Abstand zweier Punkte in erster Naherung ungedndert aft). Eq, kann als Imagindrteil 
einer komplexen Funktion aufgefaBt werden, deren Realteil dann die Spannungssumme darstellt® *. 
Spannungssumme und Drehung erfiillen die Cauchy-Riemannschen Differentialgleichungen und 
geniigen beide der Laplaceschen Differentialgleichung. Bei Platten entfallt diese Bedeutung der 
zu den Linien M, + M, = konst. orthogonalen Kurven, da wegen der fiir u und v angenommenen 
Beziehungen (4) w, identisch Null wird und auBerdem die Momentensumme an den Lastangriffs- 
stellen nicht der Laplaceschen Differentialgleichung geniigt. 

Es ist ja, wie sich aus (1), (2), (3) und (8) ableiten 1aft, 


Mz + M Ne 2 
a(n Pe) (4 gata): ie 


1 Siehe FuBnote 1 von Seite 264. 
2 Siehe FuBnote 1 von Seite 264. 
3 Siehe FuBnote 1 von Seite 260. 
4 Siehe FuBnote 1 von Seite 264, 
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Die Ausgangspunkte zur Konstruktion der Linien konstanter Momentensumme liegen am Rande, 
wo die Werte unmittelbar aus dem Isochromatenbild entnommen werden kénnen unter Beriick- 
sichtigung der jeweils giiltigen Randbedingungen. Einen Sonderfall bilden die frei drehbar auf- 
liegenden Rander. Dort verschwindet die Momentensumme langs des Randes (Naviersche Rand- 
bedingungen). a 

Aus den Randbedingungen 1a8t sich ebenfalls schon einiges iiber die Richtungen der Linien 
p = konst. am Rande aussagen. 

An einem geradlinigen, fest eingespannten Rande zum Beispiel haben die Isochromaten und 
die Linien p = konst. gleiche Richtungen, d. h. es ist dort x =f , wie folgende Uberlegung zeigt. 

An einem geradlinigen, eingespannten Rande x = konst. gilt 


Ow Ow 6? w 6?w 
athe Ses OS ee Z = ee a8 (Ne 83 
yO ney ce ae (83) 
Daraus folgt 
o 02w 
M,=—B-7=M, und M, = — Bu; =u M,=M, (84) 
und somit 
pic 85 


op Eee 1+ p aq (86) 
as, 1—p As, 
Mit (72) und (73) wird daraus 


7 cosa = 7H TP cos f. (87) 


; = und Beriicksichtigung der Abkiirzung (77) liefert 


Multiplikation mit bc ( 
a’ cosa =ccosp. (88) 
Aus der zweiten Gleichung (78) folgt dann 


Ui 


a’ cosa =e cos B =— + siny . (89) 


Dividiert man jedes Glied der ersten Gleichung (78) durch einen dieser Ausdriicke, so ergibt sich 


tga =tgh +2ctgy. (90) 
Nun ist aber y = 90°, also ctg yy = 0, und daher ist 
tga =tgB bzw. 9 "a1 = 8; (91) 


was wir zeigen wollten. 


11. Weitere Méglichkeiten zur Konstruktion der Linien konstanter Momentensumme. Wie wir 
gesehen haben, bietet die Erweiterung des Neuberschen Verfahrens die Méglichkeit, die Linien 
konstanter Momentensumme bei Platten zu konstruieren. Fiir diejenigen p-Werte, die am Platten- 
rande auftreten und die demzufolge aus den Randbedingungen und den spannungsoptischen Mef- 
werten bekannt sind, kann der Verlauf der Linien p = konst. beim Fortschreiten ins Innere der 
Platte nach (79) konstruiert werden. Die Werte auf diesen Linien sind also durch die Randwerte 
festgelegt. Zur Bestimmung der p-Werte auf denjenigen Linien p =konst., die nicht auf den 
_ Rand stoSen, dient Gleichung (80). 

Zum Abschlu8 sollen noch einige weitere Beziehungen angegeben werden, mit deren Hilfe die 

_ Ermittlung der Linien p = konst. méglich ist. 
Fallt in (35) die Richtung € mit der Richtung einer Linie p =konst. zusammen, dann ist 
Op/o§ = 0, und 6 bedeutet den Winkel, den die Linie p = konst. mit der Richtung s, des Haupt- 
biegungsmomentes M, bildet. Wir haben diesen Winkel mit ~ bezeichnet. Fiir ihn folgt aus (35) 


SAE 42 on 
tga = S2__ 92 
O40) 5 ee. (92) 
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Bequemer ist es oft, die Differentiation spannungsoptisch bestimmter Groen nach den Rich- 


tungen x, y vorzunehmen. Wir gehen deshalb von den Gleichungen (20) und (21) aus und schreiben 
sie zunachst in der Form 


op Lee wi : 8 
pit ele AGS 0) ates (q.c08 29), 
Ree (93) 


7) 7) : 
stele (q.c08 29) + 2 (qsin 2 )]. 
Diese Gleichungen setzen wir in die Beziehungen 


op Op ax ép dy op op . 
ae ve-0E by Bega gy fi8 e. 


op __ op 0x op oy _ 


Op Ox op oy Bee OR. Op 
te be Bh gg sO ap cose 


| (94) 


ein und erhalten 


Ope ea ae 
% iting 
ép_ +m 
es anya and 


7) : 7) c 
2 (qsin2g) — 2 (geos29)| + cos 0 [2 (q cos 2 @) +2 (gsin29)|| » (95) 


é - ‘ : 
px (1 8in 2 y) — © (qos 29)|—sin [2 (q cos 2 ¢) +2 (qsin2 g)||. (96) 


Diese beiden Gleichungen sind iibrigens auch zur Integration langs beliebiger Wege geeignet. 
Nun falle wieder die Richtung é mit der Richtung einer Linie p = konst. zusammen, d. h. es 
sei dp/og = 0. Der Winkel geht dann iiber in den Winkel, den die Linie p = konst. mit der 
x-Achse bildet. Wir bezeichnen ihn mit %, (es ist 0, =p +a). 
Aus (95) folgt dann 
7 0 
aa! cos 2 @) +5 (@ Bin 2) 
tg0, = aC - ——. (97) 
ae sin 2 —~) me Ce 2 9) 


Die Ausdriicke q sin 2 y und q cos 2 g werden aus den Isochromaten und Isoklinen bestimmt, 
fir hinreichend viele Schnitte x = konst. und y = konst. aufgetragen und differentiiert. Durch 
Anwendung von (97) erhalt man also zunichst ein Richtungsfeld der Linien konstanter Momenten- 
summe. Bei einiger Ubung kénnen daraus die Linien p — konst. als Kurvenziige gewonnen wer- 
den. Ganz zweckmafig wird es sein, die zu den Linien p = konst. orthogonale Kurvenschar mit- 
zuzeichnen. Die Gewinnung der Linien p = konst. aus den Richtungswinkeln gleicht etwa der 
Konstruktion der Trajektorien der Hauptbiegungsmomente aus den Isoklinen. Dies ist insbesondere 
dann der Fall, wenn zunachst die Linien oO, = konst. aufgezeichnet werden. 

Fiir solche Linien p = konst., fiir die eine Normierung durch die Randwerte nicht méglich 
ist, die also nicht auf den Rand stofen, ergibt sich eine Festlegung des Zahlenwertes aus (96). 
Diese Gleichung stellt namlich jetzt den Gradienten von p dar, den wir mit dp/ds, bezeichnet 
haben. Betrigt der Abstand zweier Linien p = konst. As,, dann folgt fiir die Zunahme Ap von 
der einen Linie bis zur nachsten aus (96): 


a 0 : rd] 
i £1 loos o> G (q sin 2y) — oy (q cos 2 9)| 


Ap. = 


z (q cos 2) + 5 (q sin 2 ®)| | As, (98) 


Am einfachsten diirfte es aber sein, zur Normierung der Linien p = konst. die Momentensumme 
durch Integration langs eines geradlinigen Weges parallel zu einer Koordinatenachse, der méglichst 
alle Linien p = konst. schneidet, zu bestimmen. Ein zweiter Integrationsweg kann zur Kontrolle 
dienen. Fiir diese Bestimmung sind die Gleichungen (22) und (23) in Verbindung mit (24) und (25) 
oder aber — was dasselbe bedeutet — die Gleichungen (93) géeignet. Aus der zweiten Gleichung (93) 
erhalt man durch Integration vom Randpunkt P, bis zam Punkte P 


—sin }, 


es ; 
1 Os 1+u 
(p)p = 5 [ a asim 20) de + (q0082 9p + Pie, 15% 2 2) dr (99) 


4 


Po 
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Die Integrationskonstante 


ia, a 
eres q cos 2 9) be (100) 


verschwindet an einem frei drehbar aufliegenden und an einem eingespannten Rande. An einem 
freien Rande wird 


2 
C=— Taye (101) 
Das eben beschriebene Verfahren ist anscheinend auch bei Scheiben zur Konstruktion der 
Linien konstanter Spannungssumme noch wenig angewendet worden. Die dazu erforderlichen 
Beziehungen lassen sich analog wie hier — ausgehend von den Gleichgewichtsbeziehungen des 
ebenen Spannungszustandes — ableiten. Man erhilt fiir den Winkel 0,, den die Linie konstanter 
Spannungssumme mit der x-Achse hildet, wieder (97), worin hierbei p die Spannungssumme 
o, +o, und q die Hauptspannungsdifferenz 0, — 0, bedeuten. Fiir den Gradienten von p ergibt 
sich bei Scheiben 


(4) ta) ta) J is 2) (7) : 
S = cos U, By (q cos 2 g) — (q sin 2 »)| + sin 0, be (q cos 2g) + es (q sin 2 9)| (102) 
mit der eben angegebenen Bedeutung von p und q. Hierbei bedeutet den Winkel zwischen der 
Richtung der Hauptspannung o, und der x-Achse. 


12. Zusammenfassung. Es sollte gezeigt werden, da sich bei der spannungsoptischen Unter- 
suchung elastischer Platten, die innerhalb der Kirchhoffschen Bedingungen auf Biegung bean- 
sprucht werden, eine ganze Reihe von Auswerteméglichkeiten ableiten lassen, die es gestatten, 
die Biegungsmomente einzeln zu bestimmen. Um spannungsoptische Bilder zu bekommen, ist es 
erforderlich, die Modellplatte aus zwei Schichten verschiedenen Materials zusammenzukleben 
(Zweischichtenverfahren) oder in der Plattenmittelflache eine reflektierende Schicht anzubringen 
(Reflexionsverfahren). Die MeBgréBen sind die Differenz der Hauptbiegungsmomente und deren 
Richtungen. So wie man beiScheibenauswertungen je nach dem vorliegenden Versuchsobjekt 
das eine oder das andere Auswerteverfahren bevorzugen wird, so hat man auch bei Platten die 
gleiche Auswahl an Auswerteverfahren. Einige davon wurden hier ausfihrlich wiedergegeben. 
Ausgehend von der Integration lings geradliniger Wege wurde die Bestimmung der Momenten- 
summe durch Integration langs beliebiger Wege beschrieben, woraus sich wiederum eine Reihe 
von speziellen Methoden herleiten lie}. Dazu gehért zunachst die Integration langs einer Haupi- 
biegungsmomentenlinie, die sich im Falle rotationssymmetrisch belasteter Kreisplatten als be- 
sonders zweckmaBig erweist, und die Integration langs einer Haupttorsionsmomentenlinie. Weiter- 
hin wurde eine Méglichkeit zur Bestimmung der Momentensumme beschrieben, die sich an ein 
von Féppl fiir Scheiben angegebenes Verfahren anlehnt, sowie ein von Neuber fiir Scheiben ent- 
wickeltes graphisches Verfahren zur Konstruktion der Linien konstanter Momentensumme bei 
Platten erweitert. Zum Schlu8 wurden noch weitere Beziehungen angegeben, mit deren Hilfe 
sich die Richtungen der Linien konstanter Momentensumme errechnen lassen. 

Der Verfasser dankt dem Leiter des Spannungsoptischen Laboratoriums, Herrn Dr. H. Schwie- 
ger, fiir fruchtbare Diskussionen. 


(Eingegangen am 22. Juni 1960.) 


Anschrift des Verfassers: Dipl.-Phys. Gotthard Haberland, Deutsche Akademie der Wissenschaften zu Berlin, 
Institut fir Angewandte Mathematik und Mechanik (Direktor Prof. Dr. K. Schréder), Spannungsoptisches 
Laboratorium, Erkner bei Berlin, FlakenstraBe 28—31. 
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Zur Inversion ebener Spannungszustande 
Von H.-J. Franeck 


1. Einleitung. Wie F. Levi Civital und J. H. Michell? gezeigt haben, bleibt die Differential- 
gleichung AAF = 0 fiir eine bestimmte Koordinatenbezichung invariant. Diese Tatsache haben 
A. Timpe® und W. Olszak* benutzt, um mit Hilfe der Inversion neue Lésungen der Bipotential- 
gleichung angeben zu kénnen. Da hierbei jedoch die Inversion in einer Abbildung innerhalb des 
kartesischen Koordinatensystems besteht, lassen die sich ergebenden Resultate unter Umstanden 
keine tibersichtliche Deutung zu. In der vorliegenden Arbeit wird daher die Inversion von Span- 
nungszustanden als Folge einer Transformation der Bipotentialgleichung in spezielle krummlinige 
Koordinaten aufgefaBt und damit eine Vereinfachung erzielt. 


2. Geometrische Grundlagen. In einer Ebene wird ein Punkt entweder durch seine rechtwinklig 
kartesischen Koordinaten x‘ oder x/ oder x* oder durch seine krummlinigen Koordinaten x* oder 
x? oder x” festgelegt. Fiir die Indizes i, j, k baw. a, fb, y ist dabei jeweils 1 und 2 zu setzen. Die 
kartesischen Koordinaten x‘ seien stetig differentiierbare Funktionen der allgemeinen Koordinaten x* 
und umgekehrt. 

Bezeichnet man eine partielle Differentiation mit einem kurzen senkrechten Strich 


EOE Gy oy eae et 


ox® oxt 
so erhalt man die Tangenten- bzw. Gradientenvektoren 
Gye i bzw. Cr 18 


und daraus bei Giiltigkeit des Summationsiibereinkommens die Komponenten des MaBtensors in 
kovarianter bzw. kontravarianter Darstellung 


ile 3 F 
Bap = 9:76 baw. gr= Od’ ci e;. 
Die kovariante Ableitung eines Tensors wird durch einen langen senkrechten Strich gekennzeichnet 


und mit dem Christoffelschen Dreizeigersymbol ['%, fiir die kontravarianten Koordinaten A* 
eines Vektors in der Form 


A*|g= AX + Av IN 6 


erhalten. Die physikalischen Koordinaten eines Tensors zweiter Stufe ergeben sich aus der Relation 


A(ep) = ee A@e, (nicht summieren) 


3. Invarianz der Bipotentialgleichung. Die Differentialgleichung 14F = 0 hat in kartesischen 
Koordinaten bekanntlich die Form 


otF ot oF . 
AAF = a(x)! + 2 a(R? a2)? % aGe)yt 0 (1) 


und geht fiir allgemeine Koordinaten in 
AAF = Fl*,°, = 6°? 6! F lapys = 0 () 
iiber. Differentiiert man (2) aus und ordnet nach den Ableitungen von F’, so ergibt sich 
AAF = F\P8y = Fjapyo 8°" 8” — Frapy 2 Dio(a"? a”? + 2 gt?) + 
+ Frage? (et? Fo. + Bi 4yo + 28% LoL e—T hol) + 
+ F.9” (P4,.(8" Doa\ie — (8 Psa)iyo] = 9- (3) 
1 F, Levi Civita, Sopra una trasformazione in sé stessa della equazione 4 a = 0, Venezia 1898 (Tip. Ferrari). 


2 J. H. Michell, Proc. Lond. M. Soc. 34 (1901/02), S. 134. 
3 A. Timpe, Z. Math. u. Phys. 52 (1905), S. 348. 
4 W. Olszak, Ing.-Archiv VI (1935), S. 402. 
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Um nun die Bedingungsgleichungen fiir die Invarianz der Bipotentialgleichung (1) zu be- 
stimmen, wird der Differentialausdruck 


> a (o F) agp BS (® P)iy411 es 2 (® F)iy122 te (@ F) e220 (4) 


ot 6—t 
angesetzt, der formal das gleiche Aussehen wie (1) hat. @ ist dabei eine Funktion der krummlinigen 
Koordinaten x1 und x2. Der Ausdruck (4) lautet dann ausgeschrieben 


Died 
2 2, Meare = (Fiaaar + 2 Farge + Fianna) ® + 

+ Fray 4 Oy + Priva 4 Dia + Finn 4 Dir + Fionn 4 Pio + 

+ Fi 2 (3 Oy, + Digs) + Fyre 8 Dire + Foe 2 (ir + 3 Diao) + 

+ Fi 4 (iar + Piee)ia + Fie 4 (ir + Piae)ja + 

+ F (Ojai + 2 Diyze2 + Pieeee) - (5) 


Aus (5) kann die Summe 
FPoyyia + 2 Pista + Foose 


berechnet und in (3) eingesetzt werden, wenn man sich bei der Wahl der allgemeinen Koordinaten 
auf solche beschrankt, die orthogonal sind (g,, = 0) und dariiber hinaus der Gleichung 


Oil oe Vg (6) 
geniigen. Man erhalt dann aus (3) 


2 2 
MAF = v= |e 2, 2, © Miasee— 
a=1 p=1 


—2 Ran( lies 2 2) Pana 7 + 2 ] 


=? Fan ( M+ 2S aes 2 Fiana( YE 2) 


Vg 
sere ay [2 (Vgis)” +2 (Vea) —Ve (Vea + Veiaa)|— 2g + 


iy mal 2 (2 + 2 (Vg) — Ve (Wei1s + Vei29)|— 2-2 


j12 Dire 1 (Dias ae Pisa Fi, 4 (Pi —— 
F erase fs 0 
@ cr 9 
und fiir die Bedingungen 
) Pi, Vgia Dye Vgie 
@ oy! 2 a i 
Vg . Vg 
310, , 0: 6 Oe Oo = = - = 
2 x Mjea 2-9 2111 x |22 __ ; 2 (Vgia)" +2 (/g\2) —V8 (V/8i1a + V/8i22)| ? “) 


Dire a Dias a Dioo0 = 0 


geht die Bipotentialgleichung fiir allgemeine Koordinaten in 
OME 
AAP = Fi*2, = — g@ a p> (D P)\ caps = 0 (8) 


liber. Die Beziehungen (7) liefern dann nur eine Lésung fiir die Funktion ®, wenn die orthogonalen 
krummlinigen Koordinaten neben der Bedingung (6) die drei Gleichungen 


5 (V/gi,)" = (/B\2)° aa /g V8i11 =0, 
5 (pe — Wea) = 4 V/s V/Sie0 = 0, (9) 
3 V8i1V8ia—2VeV E12 =0 
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identisch erfiillen. Setzt man zur Abkiirzung 

F=OF, (10) 
so ist die Bipotentialgleichung (3) fiir die den Bedingungen (6) und (9) geniigenden orthogonalen 
Koordinaten mit der Differentialgleichung 


2 * 


2 
» » Fiaxpp=0 (11) 
on al 


Li > - ° - . . i . . 
identisch. Man erhalt somit Lésungen von (3), wenn man in Spannungsfunktionen F’, die der Bi- 
potentialgleichung 


oF oF tr 


AAF = a(xt)A r 2 ACE)? (x2)? | Bt 


geniigen, die kartesischen Koordinaten x! und x? durch die krummlinigen Koordinaten x! und x2 
ersetzt. Dieser Vorgang wird als Inversion bezeichnet. Die Spannungsfunktionen 
: 


F 
F=5 (12) 


beschreiben dann neue Liésungen von Randwertproblemen. 


4. Spannungen. Die kontravarianten Komponenten des Spannungstensors werden in all- 
gemeinen Koordinaten aus der Beziehung 
Soh — e*? Bo F\y 5 
Peccinet: Mit den Bedingungen (6) und der Gleichung (12) erhalt man daraus 


Oe nay ©) QF Fise gO + F610 2 Fia(4/gDDo+ Yeo D*) + 
+ F[4)g(®,.)? —2/gOD,..— 1 OD. + /gj,0%,,]}> 
af * * = = * r 
S22 — ——__ 12 F,, /g B — Fi, (4/gOG,,4+ /g, 2) + Fels. + 
Vee re (4 /g D+ Vg. 2) - Fie \V8je (13) 


or F[4/g(®,)?—2 gO G11 + e110, Jen D,}} 
S12 — eat F4,/@2— Fi, (2/g0%,, a /Si2 2?) — Fs, (2 PD, a gi. %) Tr 
a FL gD, O,.—2/g PO y.+ /gj2P Dy + 41.0] ” 


Soll die auf einem Bogenelement ds der Koordinatenlinie x* = konst. angreifende Kraft dPé,) 
ermittelt werden, dann bestimmt man zundchst den Spannungsvektor ¢(,) in kartesischen Koor- 
dinaten. Er ergibt sich zu 


tay = Si 6; ; My) « (14) 
Fuhrt man hierin den fiir die Koordinatenlinie x* = konst. normierten Normalenvektor 
j 
nh = (nicht summieren) 
San 
und das Transformationsgesetz 
Si == oo ce 52° 
ein, so nimmt (14) die Form ops 
15 eels Say ok 
ta) = a te 


an. Mit der Scheibendicke h und dem Linienelement ds = J/g, , dx* dx? wird dann die gesuchte 
Kraft 


8a 
dPhs = > ye see SPT don : (15) 


=1 y=1 Vea 


1 4. E.Green a. W.Zerna, Theoretical Elasticity, Oxford (1954), S. 187. 
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5. Kreiskoordinaten. Die orthogonalen Kreiskoordinaten (Abb. 1) mit den geometrischen 


Werten 
1 EC x2 


_ x 
ETE oe PEEP 
1 2 x a 2 x 
t= Care etn EC 
= uae (x1)? — ()2 = = 9 xi x2 
t= d= —Pienrrayp?  S —S ee P 
L _ (ee + -_ a 
Oe eee 18 = (ap OE 


4 
11 a 


erfiillen die Gleichungen (6), (7) und (9), wenn fiir die Funktion ® die Beziehung 


Owe 


Vg 
gesetzt wird. Man erhalt dann nach (13) die 
kontravarianten Komponenten des Spannungs- 
tensors 


ye. 
Ve 
* 
sta [hat (#§—F,, x), (16) 
1 * 
shane > 7 Fie 


Die physikalischen Koordinaten werden 


* 
FE be %* * 
sq) = 42 42 (F—F,2°), 
Ig 


F Zethe® 
sea) = 7+ (FB), (17) 
§ 
- 
Abb, 1. Kreiskoordinaten. S(12) eae Pek = 
Vg 


* 
6. Beispiel. Die Spannungsfunktion F' fiir eine symmetrisch zu ihren Achsen durch konstante 
Linienlasten beanspruchte Rechteckscheibe lautet in kartesischen Koordinaten 


F=a-+ a(z)? + aa). 
Durch Inversion in die Kreiskoordinaten ergibt sich 
F = G+ a(xd)? + a(x2)? 
und damit nach Einfiihrung der neuen Konstanten a=c und a—a=c fir die Spannungen 
nach (17) 
copie 2 
S(11) = Sle—e(at], 


S22) = [e+ e@2y], ee, 
s(l2) =0. . 
Aus diesen Beziehungen ist ersichtlich, daB auf den Kreisen x1 = konst. bzw. x? = konst. konstante 


Spannungen S(11) bzw. S(22) wirken und auf diese Weise der Spannungszustand von symmetrischen, 
durch Kreise oder Kreisbégen begrenzten Scheiben unter konstanten Randbelastungen festgelegt ist. 
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ee . ° . Line 
Da fiir einen Kreis mit dem Halbmesser 7 die CréBe der allgemeinen Kreiskoordinaten durch 
a” 2 
xt = bzw. x2 = 
2r 2r 


bestimmt wird, erhalt man fiir die in Abb. 2a dargestellte Scheibe 
mit den Randbedingungen 


2 

g@=— 5: §(29)—} 
q? 

a= 5? S(22) =p 


(19) 


S(12) = 0. 


Im Berihrungspunkt der beiden Kreise wachst die Spannung S(11) iiber alle Grenzen, doch 
kann gezeigt werden, daf die in diesem Punkt tibertragene Langskraft G einen endlichen Wert hat. 
Mit den in Abb. 2b eingetragenen Resultierenden der auferen Belastung bzw. der Spannung im 

— Querschnitt x1 = 0 
> pi 


(2) 


= — 
=Yphrs Pi =2pht, pi ies 


ergibt sich die gesuchte Kraft zu 
‘ 1 2 
G=2(p—p)h 


\ 
\ 
+ 


Aus diesem Beispiel kénnen die Ergebnisse fiir den in Abb. 3a dargestellten Gelenkbogen her- 
_ geleitet werden. 
| Fir p = 0 und p= — p werden zunachst die Spannungen nach (19) 


i Sys sel r +47 (3)|> 


ma 
Pe 
pr > 2, { x7\2 
(22) =P ]1—4 # (3) 


S(12)=0. 
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Die Komponenten P(, der Resultierenden P,,) der auBeren Belastung (Abb. 3b) erhalt man 
mit den Grenzen 


=o bis x! = 00 
aus (15) zu Z 
af R2} es Rr 
Poa Oph Fase Me Oe cate 


Abb. 3. Gelenkbogen. 


Die Resultierende hat die GréBe 


1 


R 
Be pee 
/ R? 4} 
Analog ergeben sich mit den Grenzen 
1 a ; 2 a* 
os ~ bis ox? = 
bg 27 E 


fiir die Komponenten Hee der resultierenden Stiitzkraft P.,, 


= rm rfr—R 3 Rr? 
r i Be x 
Pepa Shea Foe P? = 2 ph——_,- 


re R?+ 22 
Poa teh yy Vt 
Damit folgt fiir die Gelenkkraft G aus Gleichgewichtsgriinden 


1 2 


Ca aph— =: 
P r+? 


7. Zusammenfassung. Die homogene Bipotentialgleichung fiir eine Spannungsfunktion F 
wird in allgemeinen Koordinaten angeschrieben. Es l4Bt sich zeigen, daB diese Gleichung unter 
ganz bestimmten Bedingungen formal das gleiche Aussehen wie die Bipotentialgleichung in ortho- 
gonalen kartesischen Koordinaten hat. In diesem Falle kénnen aus bekannten Lésungen in karte- 
sischen Koordinaten neue Ergebnisse durch Inversion erhalten werden. 


(Eingegangen am 1. August 1960.) 
Anschrift des Verfassers: Dr. rer. nat. Heinz-Joachim Franeck, Dresden A 27, Am Gericht 20. 
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Die Scheibe mit elliptischem Kern 
Von G. Kaiser 


1. Einleitung. Eine unendlich ausgedehnte Scheibe konstanter Wandstarke enthalte einen durch 
eine Ellipse begrenzten Bereich, der mit einem Material eines anderen Elastizitatsmoduls ausgefiillt 
sei. An der Grenzkurve seien beide Stoffe zug-, druck- und schubfest miteinander verbunden. 

In der vorliegenden Arbeit! wird der Spannungszustand beschrieben, den diese Steifigkeits- 
anderung hervorruft, wenn die Scheibe einen homogenen Grundspannungszustand unterworfen 
wird oder wenn sie — bei unterschiedlichen Warmedehnungskoeffizienten von Innen- und AuBen- 
raum — einer Temperaturanderung ausgesetzt ist. Fiir einige einfache Grundspannungszustande 
werden die Maximalspannungen in Diagrammform angegeben. Dehnungsmessungen an einer 
Scheibe mit einer elliptischen ebenen Verstarkung bestatigen die Brauchbarkeit dieser Ergebnisse 
fir die in der Technik haufiger anzutreffenden drtlichen Verstarkungen, bei denen also der Steifig- 
keitsunterschied durch die verschiedene Wandstarke hervorgerufen wird. 


2. Spannungs-Verzerrungs-Beziehungen in elliptischen Koordinaten. Das elliptische Koordi- 
-natensystem hangt mit dem cartesischen durch 


x =ccoshacosf, (la) 
y = csinha sin 6 : (1b) 
zusammen. Dabei ist 2c der Brennpunktabstand, und durch « = konst. sind die Ellipsen 


x2 yy 


c? cosh? « a sinh? a Ee 
durch B = konst. die Hyperbeln 
Be gu eS 1 
c? cos? B c? sin? B 


bestimmt. Der Faktor h der Verzerrung zwischen den Elementen des cartesischen, dx und dy, und 
denen des elliptischen Koordinatensystems, da und df, lautet? 


2 
2 2 
oe 6 =e (2) 


Zwischen den Verschiebungen u, und ug, den Dehnungen é,, ég und é,,, der ebenen Dehnungs- 
summe e’ = €, + é, und der Drehung w bestehen mit 


We __ ug 

en. und Vier 

die Beziehungen 
ou u Gh? v oh? oa 
eta Ton 2 aR ie) 
év v oh? u oh? 3b 
eg 0 BRS a ee) 
(7) (a) 

ean = £0) + Gy (Pw), (30) 
aN DRE de 3d 
at e + 6B)? oy) 
_ Hh (ov ou 3e) 


i at fii : i hschule Hannover geneh- 
i! iirzte F einer von der Fakultat fiir Bauwesen der Technischen Hoc 
eee isseatatio cites: Prof. Dr.-Ing. habil. A. Pfliiger, Korreferent: Prof. Dr. rer. nat. H. Kauderer. 
2 C. E. Inglis, Transactions Institution Naval Architects 55 (1913), S. 219. 
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Aufer den Spannungs-Dehnungs-Beziehungen des ebenen Spannungszustandes 


og = 26 (e+ Ere (4a) 
lt j 
lig SG (40) 


in denen G den Schubmodul und y die Querzahl bedeutet, werden noch die Gleichgewichtsbedin- 
gungen bendtigt, die hier in der durch Dehnungssumme e’ und Drehung @ ausgedriickten Form 


angegeben werden sollen: 


de’ ow 

ba = (1 7) ee 0, . (5a) 
de’ 0 

gt te) = 0. (5b) 


Diese Cauchy—Riemannschen Differentialgleichungen sind bekanntlich von jeder Funktion der kom- 
plexen Veranderlichen a + if erfiillt, wenn man ihren Realteil als e’ und ihren Imaginarteil als 
(1 —) w auffaBt. Ausgehend von zweckmaBigen Ansatzen fiir solche Funktionen lassen sich die 
Differentialgleichungen (3) integrieren, so daB man Verschiebungen und Spannungen in Abhangig- 
keit von spater durch Randbedingungen zu bestimmenden Freiwerten erhalt. 


3. Zu beiden Achsen symmetrische Spannungszustaénde. Innerhalb der Grenzellipse mu8 die 
Funktion fiir « = § = 0 (in den Brennpunkten) stetig mit stetigen Ableitungen sein. Der Ansatz 


F(« +16) = Ccosh n(a +76), (6) 
worin C eine Konstante und n eine ganze Zahl ist, fiihrt tiber Real- und Imaginarteil 

Re F= e’ = Ccoshnacosnf, (7a) 

Im F = (l—yp) » = Csinh na sinnf (7b) 


und die Gleichungen (2) und (3) auf das Differentialgleichungssystem 


i ig = CE {[eosh (n + 2) « + cosh (n — 2) a] cosnf 
— cosh n « [cos (n + 2) 8 + os (n — 2) B]}, (8a) 
wast 5 =5 (= {[sinh (n + 2) « + sinh (n — 2) a] sinn B 
— sinh no [sin (n + 2) 6 + sin (n — 2) ]}, (8b) 
dessen allgemeine Lésung mit den Abkiirzungen 
_3—4p anes Gre 
ay < r und A, => 16 (1 — jp) 


und der beliebigen Konstanten B,, lautet 


AW ee , . 
u= dn |(7 7 1) sinh (+ 2) « cos nea 1) sinh (n — 2) « cos np 


n+ 1 
A ; A 
~(—5 se 1 sinh no cos (n + Ee eats 1) sinh n a cos (n—2) 6) 
+ B, sinh na cosnB, | (9a) 
ree 5 
PS 4n (544+ 1) cosh (n + 2) «4 sinn 6 + a 1) cosh (n — 2) a sinn B 
A : A 
ae eae cosh na sin (n + Jeers 1) cosh n a sin (n—2) | 
— B, coshnasinnf. (9b) 
Einsetzen in (3) und (4) liefert nach einiger Zwischenrechnung mit den Abkiirzungen 
26 4, 26 n, 
An = ~= An, B, => B, 
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die Spannungen 
Ox (cosh 2 « — cos 2 8)? = 
= A, {(2 + n) cosh na cos (n + 4) 6 
— [4 cosh (n + 2) a+ (4+ 2 n) cosh (n — 2) a] cos (n + 2) B 
++ [(2 —n) cosh (n + 4) 0 + 8 cosh na + (2 + n) cosh (n — 4) a] cos nf 
— [(4 — 2 n) cosh (n + 2) a + 4 cosh (n — 2) a] cos (n — 2) B 
++ (2 —n) cosh n a cos (n — 4) B} 
+ B, {— (n — 1) cosh n a cos (n + 2) 6 
++ [(2 — 1) cosh (n + 2) « + (n + 1) cosh (n — 2) a] cos np 
— (n + 1) cosh na cos (n — 2) 8}, (10a) 
6g (cosh 2 « — cos 2 f)? = 
= A, {(2 — n) cosh na cos (n + 4) B 
— [4 cosh (n + 2) a + (4 — 2 n) cosh (n — 2) a] cos (n + 2) 
+ [(2 + n) cosh (n + 4) « + 8 cosh na + (2 — n) cosh (n — 4) a] cosnB 
— [(4 + 2 n) cosh (n + 2) « + 4 cosh (n — 2) a] cos (n — 2) B 
+ (2 + n) cosh na cos (n — 4) B} 
— B, {— (n—1) cosh n a cos (n + 2) B 
+ [(n — 1) cosh (n + 2) « + (n + 1) cosh (n — 2) a] cos nf 
— (n + 1) cosh n a cos (n — 2) 8}, (10b) 
Txg (cosh 2 « — cos 2 f)? = 
= A, {— nsinh na sin (n + 4) 6 + 2 n sinh (n — 2) a sin (n + 2) 8 
+ n [sinh (n + 4) a — sinh (n — 4) a] sinn B 
— 2 nsinh (n + 2) a sin (n — 2) 6 + n sinh na sin (n — 4) B} 
+ B, {— (n— 1) sinh na sin (n + 2) B 
+ [(n — 1] sinh (n + 2) a + (n + 1) sinh (n — 2) a] sinn B 
— (n+ 1) sinh na sin (n — 2) B}. (10c) 


Mit geraden n und beliebigen 4, und B, stellen die Gleichungen (9) und (10) mégliche Verschiebungs- 
und Spannungszustande in einer Scheibe mit auBerer elliptischer Begrenzung dar. 
Fir eine unendlich ausgedehnte Scheibe, die im Innern durch eine Ellipse begrenzt ist, also fiir 
eine Scheibe mit elliptischem Loch, hat Inglis* Spannungs- und Verformungszustande angegeben. 
Aus diesen unendlich vielen Lésungen werden fiir den Innenraum (Index 1) die fiir n = 0 und 
n= + 2, fir den AuBenraum (Index a) die fiir n = + 1, —1 und —3 ausgewahlt und addiert. 
Im Kern lauten diese Anteile der bezogenen Verschiebungen und Spannungen 


26; 


= u; = B, cos 2 6 sinh 24 + 2(A—1) Ay sinh2a, (11a) 
26; 
a v, = sin 2 8 [—2 (A— 1) A, — B, cosh 2a], (11b) 


Oy,; (cosh 2 « — cos 2 f)* 
= cos 4 6 [4 A, — B, cosh 2 a] 
+ cos 2 8 [— 16 A, cosh 2 « + (3 + cosh 4 a) B,] 
+ (8+ 4 cosh 4a) A, — 3 B, cosh2a, (11c) 
6,,; (cosh 2 a — cos 2 f)? 
= cos 4 [4 A) + B, cosh 2 a] 
+ cos 2 B [— 16 Ay cosh 2 « — (3 + cosh 4 «) B,] 


+ (8 + 4 cosh 4a) Ay + 3 B, cosh2a, (11d) 
Txp,¢ (cosh 2 « — cos 2 f)? 
= B,sin4f sinh 20 + sin 2 6 [— B, sinh 4a], - (116) 


* Siehe FuGBnote 2 von Seite 275. 
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Es zeigt sich, daB sich mit den Gleichungen (11) und den von I nglis angegebenen Liésungen des 
AuBenraumes alle Forderungen fiir die Grundspannungszustande des einachsigen Zuges parallel 
zur kleinen oder groBen Halbachse befriedigen lassen, sowohl die Randbedingungen im Unendlichen 
als auch die Ubergangsbedingungen an der Grenzellipse zwischen Kern und Aufenraum. An der 
Grenzellipse « = 4 = konst. mu8 innen und aufen Gleichheit der Verschiebungen gefordert wer- 
den, da nach Voraussetzung keine Klaffungen auftreten kénnen, und Gleichheit der Normalspan- 
nungen im elliptischen Schnitt und der Schubspannungen, da Gleichgewicht am Element vorhanden 
sein muB. Diese vier Bedingungen 


(uj)z = (Ua) » (12a) 
(v:)5 = (Ya)a » (12b) 
(64,1) = (a,a)a » (12c) 
(Tapia = (Tap,a)a (12d) 


fiihren, da sie fir § identisch erfiillt sein miissen, auf acht Bedingungsgleichungen: 


a 
Ay Bis By ee Big Rechte Seite 
EEE EEE IEEE SEER 
—4 cosh2a% | 0 2 | e2% —17 
16cosh20%  |—3—cosh 4a 0 one cosh 2 & bee (3 + e-44) 4 cosh 2a + e¢* 
—8—A4cosh4«| 3cosh2« —2 sinh. 2% 4+ 26-40 | 3 2% | —2—et® + 3 e248 
0 — sinh 2 0 | 0 ea | + 2% 
0 sith 41 2 7) desk a) GE 2) | = 2 ah oo 
0 — sinh 2a 0 | e(1 +A) eee | Se(l +A) + eet 
=o 1} | ‘ = = 
sinh 2 % 0 +e e (1 —A) e 2% 0 gel — A) e2% 
2(a—1) cosh 2 & 0 e(1 —A) ee 2% — Fela # eet 
(13a—h) 
Dabei ist ¢ = G,/G, das Verhaltnis der Steifigkeiten; das mena, Vorzeichen gilt fir Zug 


parallel zur aaa Halbachse. 


Da nicht alle Gleichungen voneinander unabhangig sind, geniigen die fiinf Konstanten zur Er- 
fillung der acht Bedingungen; weitere zwei Konstanten folgen aus der Erfiillung der Randbedin- 
gungen im Unendlichen, wo die Spannungen dem Grundspannungszustand anzupassen sind. Die 
sieben Konstanten haben die Werte 


Ay=H {ze + Illes + 1) e#—(—e4¥] Fo(e—1) A+ I, (14a) 
B, =p 26+ 1) [I G—-1) FRe+i— ea], (14b) 
Be =y—1) (A— 1) {— (eA +1) et +e (A+ 1)4 (€—]) [e444 4sinh2a]}, (14e) 
An =ylreGt 1)? e2@ + (¢—1) Q ée+A—1) (1+ 2 5) sinh 23 : (14d) 
By = FF {e(e—1) (A+ 1) (A—1) 2 sinh 2a 
: + [e (4 + 1)? — (e— 1) (2e +A—1) €24 2 sinh 2 a]}, (14e) 
und 
Bj= Fl (14g) 


mit dem Nenner N = (€4 + 1) (2¢ + 4—1) 8% — (ec — 1) (2eEA—A+ 1) e232, 
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Der Spannungszustand, ausgedriickt als Vielfaches des Grundspannungszustandes des einachsigen 
_Zuges wird durch die folgenden Gleichungen vollstandig beschrieben: 


Innenraum: 
Ox B 3 h 4 2 6 — cosh 2 4 
aU erate 2 ____ cosh 4 ¢ cos 2 6 — cosh 2 « cos 46 
OS 9 8 cosh 2 « — cos 2B (cosh 2 « — cos 2 f)? ? (15a) 
op On, 
ee ae (15) 
Top B, sinh2¢sin2 8 
Cope 4 cosh 2«—cos2f ’ (15¢) 
AuBenraum: 
On I 
Gx. 8 (cosh 2 « — cos 2 f)? 
{cos 48 [—2 A_,—e*B,4+2 Aaa es? eBe| 
+ cos 2 8 [8 A_, cosh 2a + (3 + e#%) B.,—8 A,, cosh 2.0 
sel Rae OS SES 
— (4+ 2 e+) A ,—3 eo B,— 2 B_, sinh2a 
oe (Ae eA 3 ee Bt, (15d) 
CRs 1 
Oxo 8 (cosh 2 @ — cos 2 f)? 
{cos4B[—2A_,+ 2*B,+ 2 Aj —e¢ 25 B.,| 
+ cos 2 8 [8 2% A_, — (3 + e44) B_s— 8 e-24 A,, + (3 + e—42) B,,] 
— (4+ 2 e+) A ,+3e«B 4+ 2 B_,sinh2« 
+ (4 +2 e 4%) A,—3e 2+ BB}, (15e) 
Ta B ioe 1 
Oo. 8 (cosh 2 « — cos 2 f)? 
{sin 4 6 [e?* B_,-+ e—2o Bis] 
+ sin 2 8 [— 4 A_, cosh 2a — (3 + e4*) B ,—2B, 
— 4 A,, cosh 2 « — (3 + e—44) B,J}. (15 f) 


In 4hnlicher Weise laBt sich der Spannungszustand in einer Scheibe angeben, deren elliptischer 
Kern aus einem Material mit anderer Warmedehnungszahl «7 besteht und die einer Temperatur- 
anderung AT ausgesetzt ist. Die rotationssymmetrische Temperaturverformung einer Scheibe 
lautet in elliptischen Koordinaten, ausgedriickt durch die bezogenen Verschiebungen, 


up = Sap AT sinh 2a, (16a) 
vp = —S ap AT sin 26. (16b) 


Im Innenraum wird der Stérspannungszustand durch die Gleichungen (11) beschrieben. Aus den 
Randbedingungen im Unendlichen, o, = og = txg = 0, und den Ubergangsbedingungen 


U; + Up; = Ua + UT, > (17a) 
v; + UT,; = Ye + YT a> (17b) 
On,i a Oa,a 2 (17 c) 
3 Tx B,i a Tx B,a (17d) 
lassen sich die Konstanten errechnen: 
Ay= xq lle + 1) oF + (0-1 eI], (18a) 
B,=y2et+1), (18b) 
PA (18c) 


B,=0, 3 (18d) 
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B,=—p[leA + 1) ##—e(Q+ 1) + @—Ne*a], (186) 
Ay = — py (€—1) (#* —e*8), (18f) 
Bi = we (A+ 1) &% (@% — 88), (18g) 


Zur Beschreibung des Stérspannungszustandes kénnen die Gleichungen (15) zusammen mit den 
Konstanten nach (18) benutzt werden, wenn man dort o,, durch die BezugsgréBe 


or = G6; A4T (a7, 4 aa Xr, i) (19) 


ersetzt. 


4, Zu den Achsen antisymmetrischer Spannungszustand. Fiir reinen Schub wablt man die i-fachen 
komplexen Ansatzfunktionen und gelangt auf einem dem Rechengang des vorigen Abschnitts ent- 
sprechenden Wege zu den Spannungen im Innenraum 


o. _ D, sinh2asin26 (20a) 
Tzy,c 2 cosh2a—cos2p’ 
geen (20b) 
EN pen eet BP — 4 ita eT eee eee (20¢) 
Tzy,0 4 cosh2a—cos2f (cosh 2 « — cos 2 f) 


und im AuBenraum 


Caen 1 
Try,  4(cosh 2 ~ — cos 2 f)? 


X {sin 4B [2 C, + e—24 D, — eo] 
+ sin 2 6 [— 2 C; cosh 2 a — (3 + e—4*) D, + 3+ ef4]}, (20d) 
i} Soe 1 
Try, 4 (cosh 2 «— cos 28)? 
X {sin 4 6 [2 C, — e—2« D, + e2«] + sin 2 6 [(3 + e—4%) D, —3 — et4]} , (20e) 


ih es 1 
Try.~  4(cosh 2 « — cos 2 8)? 


X {cos 4B [— e—24 D, — e24] 


+ cos 2 8 [C, sinh 2.0 + (3 + e—44) D, + 3 + eta] — 3 e—24 D, — 3 e2 a} (20f) 
mit den Konstanten 
C= ee eee = 2la 
2 GAL VD) 8? es ete Gre) 
2e(A + 1) e2* (21b) 
(eA + 1).e?* + (ce — I) e 24 

> (g— 1) 628(e2% _ 5-28 
pee Aero eee (21c) 

(eA + 1) e?* + (e— le 24 
Dos e*® [(¢ — 1) e&* + (eA + 1) e 25) (214) 


(cA + Ve?* + (— 1) 62% 


5. Auswertung. Mit den in den Abschnitten 3 und 4 aufgestellten Beziehungen ist es méglich 
den Stérspannungszustand einer Scheibe zu berechnen, die im Unendlichen einem homogenen 
Spannungszustand unterworfen ist oder, anders formuliert, die bei gleichen Elastizitatseigenschaften 
von Innen- und AuBenraum einem beliebigen homogenen Spannungszustand unterworfen ware. 
Man hat dafiir den Grundspannungszustand in seine Komponenten o,.,., ,,, und T, yoo Zu zerlegen, 
darauf die Gleichungen (14), (15) und (20), (21) anzuwenden und die so berechneten Spannungen 
zu summieren. Mit A = (3 —u)/(1 + y) gelten diese Formeln fiir den ebenen Spannungszustand; 
fiir den ebenen Verzerrungszustand sind die gleichen Beziehungen giiltig, jedoch ist A = 3 — 4 be 
zu setzen, 
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Aus den Gleichungen (15a, b, c) laBt sich nachweisen, da der Spannungszustand im Innenraum 
fiir einachsigen Zug homogen ist. Der mit A, behaftete Lésungsanteil stellt den Zustand o, = Oys 
Txy = 0, der mit B, behaftete den Zustand o, = — Oy, Tey = 0 dar. Auch unter reinem Schub 
stellt sich ein homogener Zustand, namlich wieder reiner Schub, im Innenraum ein, wie sich aus 
(20a, b, c) herleiten laBt. Es ist damit das bemerkenswerte Ergebnis gewonnen, da8 der Spannungs- 
zustand im elliptischen Innenraum immer homogen ist. Auch unter einachsigem, beliebig gerichteten 
Zug entsteht im elliptischen Kern ein vom Ort unabhangiger Spannungszustand, der allerdings im 
allgemeinen nicht mehr einachsig ist. Das Abklingen oder Anwachsen der Spannungen beschrankt 
sich also nur auf den AuSenraum. 

Erstaunlich ist, daB diese Tatsache, die natiirlich auch fiir den Sonderfall des Kreises gilt, bei 
der Berechnung des kreisférmigen Kernes von Goodier! unerwahnt geblieben ist. Erst bei Muskhe- 
lishvili? und M. Davin® finden sich Hinweise darauf. Diese Homogenitat ist auch nicht etwa eine 
Folge der hier getroffenen Annahme gleicher Querdehnungszahlen fiir Innen- und AuSenraum, son- 
dern gilt allgemein fiir homogene Grundspannungszustande. 


Abb. 2. Bezogene Hauptspannung 0,9, infolge ox oo: 


ier, J. appl. mech. 1 (1933), S. 39. ; © : 
: Nn 7 Sracthclishen, Sone bed pine of the mathematical theory of elasticity, Groningen 1953. 
5 M. 


if 
Davin, Archiwum Mechaniki Stosowanej 2, XI, Warschau 1959, 
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f=43 
07=G,4 (cer 9 Gr) 
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infolge Temperaturdehnung. 


Abb. 3. Bezogene Hauptspannung o 


max 


Abb. 4. Bezogene Hauptspannung Cmax infolge Spoor P = 30°. 


Abb. 5, Bezogene Hauptspannung o,,,,, infolge a, 


pog? P= 7730", 
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In Abb. 1 bis 5 ist fiir mehrere Grundspannungszustande das Verhaltnis der gréBeren Haupt- 
spannung zur Bezugsspannung in Abhangigkeit vom Ort aufgetragen. Es ist stets das Halbachsen- 
verhaltnis a/b = 5, das Steifigkeitsverhaltnis ¢ — G,/CG, = E,/E, = 3 und uw = 0,3 zugrundegelegt. 
Bei den drei doppelsymmetrischen Spannungszustanden allseitiger Zug, Temperaturdehnung und 
einachsiger Zug parallel zur groBen Halbachse ist die im gesamten Kern herrschende Spannung 
gleich dem GréBtwert der Hauptspannungen des AuBenraumes. Dagegen liegt fiir einachsigen Zug 
in Richtung von 30° zur groBen Halbachse die Maximalspannung nicht mehr auf der groBen Halb- 
achse und ist geringfiigig gréBer als die Spannung des Kernes. 


max Sax 
50 7 == = 
30 
20 4 | r | 
10 
[ 
8 Zp af 
y eZ 

6 KD) 

40 
4 +f 15 
3 5 

im 

3 
é 

é 

ii 

& 

4 6b 
7 é 3 4 § 6 789” 20 4Q 60 100 


Abb. 6. Bezogene GréStspannungen infolge allseitigen Zuges in Abhingigkeit von a/b und «. u = 0,3. 


Fiir die doppelsymmetrischen Spannungszustande kann fiir ¢ > 1 der GréBtwert der gréBeren 
Hauptspannung leicht in Abhangigkeit der Halbachsen- und Steifigkeitsverhdltnisse angegeben 
werden, indem man in (15a) 6 = 0, « = a und e24 = (a + b)/(a — b) einsetzt, namlich 
fir allseitigen Zug 


b —b 
: (4+) 2424 @€—NG—)—C—NS 
fm e+) oe Se) 
CA @et4—1)- +> —(@—}) @si—14 1a 
fir Zug parallel zur groBen Halbachse 
b —b 
ee ee 9 Ge me 
a b —b 
Or00 (A+ I) Re t4—) LP —(—N@ea—A + NA 
und fir Temperaturdehnungen 
a+b a—b 
ams 
eel (A+ 1) teAt)t+€ es (220) 
a+b a—b 


FANG +4—YSt_@—N@ea—2+ 0° 


a—b +5 


Sie sind in Diagrammform in Abb. 6, 7 und 8 aufgetragen. 
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Fiir einige Grenzfalle ergeben sich recht einfache Ausdriicke, die hier noch angegeben seien. 
Bei ¢ = 00, einem starren elliptischen Kern, lauten die Achsenrandspannungen 0, fiir allseitigen 


Zug o,,, Spannung auf ay Halbachse 


kleiner 
Gaal A Lames a/b 2 
Go «4A ce Call ie 


Abb. 7. Bezogene GroStspannungen infolge o, ,, in Abhingigkeit von a/b unde. uu = 0,3. 


2 roBen B 
fir Zug o,, parallel zur ae Halbachse, Spannung auf ee Halbachse 
Ox +) b 
ee ECE |. (23b) 
fii kleinen groBer 
ir Zug o,, parallel zur te Halbachse, Spannung auf eee Halbachse 
On os A + i a/b 
a= aE + (A— 3) ale (23) 


Die von Inglis bereits berechneten Randspannungen auf der grofBen Halbachse eines elliptischen 
Loches folgen aus (15b) durch Einsatzen von ¢ = 0,6 = 0 und a = a. Sie lauten fiir allseitigen Zug 


C3 ows, a ; 
— = 2- 5? (24a) 
fiir Zug parallel zur groBen Halbachse 


Sinan dl (24b) 
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fiir Zug parallel zur kleinen Halbachse 


(24c) 


Labt man a/b iiber alle Grenzen wachsen, so kommt man zu einem unendlich schmalen, einem 

Bee igen Kern. Der Wert b = 0 liefert, in die Gleichungen (22a—c) eingesetzt, fiir allseitigen 
ug 

; Ox hal 3—A 


=) € | 
Ope ae ay ea (25a) 
GO, 
hess Ses 
04 aoe 
83 
92) ah 
19 
of 
an 
Q1 i 
Ay 
49 
005) AY zi ue 
aif — 
5 rie | 
Ou Q le aS 
Ay 
001 
io + 
2 bs cele ee 
0,005) 
= 
if 
Ay Z| 
0,002 | + 
a 
B 
G00, Cat Gea S68 8 20 1 60 700 
Abb. 8. Bezogene GréStspannungen infolge Temperaturdehnung in Abhingigkeit von a/b unde. u = 0,3. 
fiir Zug parallel zur groBen Halbachse 
; to 
—“=g¢ (25b) 
Too 
und fiir Temperaturdehnungen 
Ou e—]1 1 ( 
25¢) 


oe oar Iv? eel 
Die Gré8tspannungen fiir den Sonderfall des kreisférmigen Kernes betragen fir allseitigen Zug 


Gar e(A+1) 
ey Y en) 
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fiir einachsigen Zug 


Ox _ e(A+1)(A+2¢ +4) (26b) 
Oo 2(eA+1)(2e+A—]) 
und fiir Temperaturdehnungen 


La aero | (26¢) 
Gp. 2s k= 1) 


6. Dehnungsmessungen an einer elliptischen Verstarkung. In technischen Anwendungen treten 
Versteifungen 6fter in Form laschenartiger Verstarkungen auf. AuBer den dabei entstehenden 
Kerbspannungen am Wandstarkensprung bildet sich ein Stérspannungszustand aus, der sich fir 
nicht zu groBe Unterschiede der Wandstarken t; und t, gut mit den entwickelten Gesetzen beschrei- 
ben laBt. Das Steifigkeitsverhaltnis ist dann ¢ = G;#,/G, t,, und die Gleichungen liefern statt der 
Spannungen die Langskrdafte. 


= 
Ee 
SEES CRETE z 
eau —_—_—_S_ 
aw _ 
| es 


0 100 200mm 


Abb. 9. Spannungserhéhung an einer elliptischen Verstarkung. 


Durch Dehnungsmessungen an einem Versuchsstiick aus Leichtmetallblech, auf das zwei ellip- 
tische Bleche gleicher Wandstarke aufgeschweifSt waren, wurde diese Annahme bestatigt. Das 
Halbachsenverhiltnis, bis zur SchweiBnahtmitte gemessen, betrug a/b = 4,1. Nach Abbau der 
Eigenspannungsspitzen durch Uberbelastung warden die Dehnungen unter einachsigem Zug parallel 
zur groBen Halbachse gemessen. Einen Vergleich der MeSpunkte mit der Rechnung zeigt Abb. 9, 
in der der Verlauf der Spannung o, langs der x-Achse dargestellt ist. Die Ubereinstimmung im AuBen- 
raum ist recht gut, im Innenraum ist die auf der Oberflache der Verstarkungen gemessene Dehnung 
jedoch etwa 7% geringer als die errechnete. Die Rechenvoraussetzung der konstanten Spannung 
tber die Wandstarke trifft infolge der Elastizitat der Schwei8naht und der raumlichen Spannungen 
hier nicht zu. 


_ 7. SchluBbemerkung. Fiir belicbige homogene Grundspannungszustande werden die Stérspan- 
nungszustande berechnet in Scheiben, die einen durch eine Ellipse begrenzten Kern einer anderen 
Dehnsteifigkeit enthalten. Von dieser Lésung ausgehend la8t sich auch der Spannungszustand 
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angeben, der durch einen Kern mit unterschiedlichem Warmedehnungskoeffizienten bei Temperatur- 
anderung entsteht. 


Es ergibt sich, da unter allen homogenen Grundspannungszustanden der Spannungszustand 
im Kern immer homogen ist. 


Aus einigen Diagrammen lassen sich die an einem elliptischen Kern auftretenden GréBtspan- 


nungen fiir drei einfache Grundspannungszustande in Abhangigkeit von Halbachsen- und Steifig- 
keitsverhaltnis ablesen. 


SchlieBlich wurde noch gezeigt, daB sich die Ergebnisse auf die in der Technik anzutreffende 
értliche Verstarkung einer Scheibe gut iibertragen lassen. 


(Eingegangen am 8. August 1960.) 


Anschrift des Verfassers: Dipl.-Ing. Gottfried Kaiser, Miilheim/Ruhr, Saarnberg 71. 
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Die Schwingungszahlen linear-elastischer Systeme 


Von O. Bottema 


1. Vor kurzem hat Neuber! bei linear-elastischen Punktmassensystemen gewisse GesetzmaBig- 
keiten fiir die Frequenzen gefunden und auf die Berechnung von Frequenzschranken angewandt. 
Ein solches System besteht aus n Massenpunkte A; (i = 1,..., 1) im Raum, wo 4; die Masse m; 
hat und A; und A; (i ~ j) durch eine lineare Feder mit der Federkonstante ¢;; (= c;;) verbunden 
sind. Der grundlegende Gedanke der Untersuchung bestand darin, daB anstelle der Verschiebungen 
die inneren Krafte als Unbekannte eingefiihrt werden. ,,Die zur Errechnung der Eigenfrequenzen 
dienende Determinante wird so wesentlich einfacher und ibersichtlicher“. Die Federkonstanten 
treten nur in den Diagonalfeldern auf, und dieser einfache Aufbau gestattet, allgemeine Gesetze 
fiir die Frequenzen aufzustellen. 

Wir zeigen in der vorliegenden Notiz, daB, auch wenn man die Verschiebungen als Unbekannte 
beibehalt, eine Frequenzdeterminante entsteht, woraus man in einfacher Weise die erste Neubersche 
Formel ableiten kann. Daneben bemerken wir, daB diese Frequenzdeterminante im allgemeinen 
einen Vorzug hat vor derjenigen, die aus einer Betrachtung der inneren Krafte entsteht. Ist das 
System eindimensional, so hat es n Freiheitsgrade, ein zweidimensionales hat 27, ein dreidimensio- 
nales 3n Koordinaten. Daraus geht hervor, daB die Anzahl der Frequenzen in diesen drei Fallen 
baw. gleich n— 1, 2n—3 und 3 n— 6 ist, da ja die Bewegung des Systems, als ein Ganzes be- 
trachtet, die Freiheitsgrade eines starren Kérpers aufweist. Die Frequenzdeterminante, welche aus 
der Betrachtung der Verschiebungen entsteht, ist dabei von der Ordnung n, bzw. 2n und 3n (und 
die beziigliche Gleichung hat eine einfache bzw. dreifache und sechsfache Wurzel Null). Geht man 
aber von den inneren Kraften aus, so entsteht eine Frequenzdeterminante, deren Ordnung der 
Anzahl dieser Krafte gleich ist; diese Anzahl betragt aber, wenn samtliche c;; 4 0 sind, nicht 
weniger als 1/2 n (n — 1). Wir zeigen weiterhin an einem Beispiel, daB die zweite Neubersche Formel 
in der angegebenen Fassung nicht immer zutrifft. 


2. Wir betrachten vorlaufig den Fall eines ebenen linear-elastischen Systems, gehen von einer 
Gleichgewichtslage A? aus, wobei samtliche Federn ihre natiirliche Lange haben, und wahlen in der 
Ebene ein rechtwinkliges Koordinatenkreuz Oxy. Der Winkel des Vektors 4?A? mit der x-Achse 
sei yj; man hat cos y;; = — cos g;;, sin ~;; = — sin p;;. Die kleinen Verschiebungen des Massen- 
punktes A; aus seiner Gleichgewichtslage seien in der x- und y-Richtung bzw. x; und y;. Dann 
gilt fiir die kinetische Energie T und die potentielle Energie V des Systems 


1 : 4 
T=> Sim (i+), 


1 . 
i eran [(%; — x;) cos Piz + (yj— y;) sin Vij)?» | 


i>y 


(1) 


und die Lagrangeschen Gleichungen fiir die Koordinaten x, und y, sind also 


™y, %, + » in [(x, — %;) cos* pi, + (¥, — ¥;) cos Y;, sin py] = 0, 


; 2 
m, Yy_ DH Cin [(%_ — %;) Cos v4, sin p;, + (y;,— ¥;) sin® Pil = 0. ee 


Die Frequenzgleichung ist somit 


D G1 608" Yjy — mM, wD) c;1 Cos H;, SiN Y; — €y1 COS? Poy — €g1 COS Hg SIN Yo, 
| 2% COS G1 SiN G;,  D' ¢;, sin? y;; — Mm, w? — cg, COS Y_, SiN Po, «2s ese waa 

0= | 4s COS" 5 = %_ COS Py 2SIN Pg D Cig COS? Y;.—mMygw* .....45..... -(3) 
—— &12 608 Pye SIN P15 a Gyg Bt Og SiS ere Gece ne ete eee ae a 


1 H. Neuber, Ing.-Arch. 28 (1959) S. 229. 
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und wir erhalten eine Gleichung 2 n-ten Grades fiir w? mit den Wurzeln «?. Teilt man noch die erste 
und zweite Zeile der Determinante durch m,, die zwei folgenden durch m, usw., dann zeigt die 
Betrachtung der Glieder der Hauptdiagonale: 


2n 
2 oh = 2, Bee 65 008% gj + DY) cj sin? si sa eS it 7) e 
= J 4 : “i 


Nun ist aber 
= iloarpaat 
Vien ay eS a3 i) (5) 


das Frequenzquadrat des Zweimassensystems A;A;, das durch Herauslésen aller Federn bis auf 


das Verbindungselement zwischen A, und A; entsteht. Damit ist das eine Neubersche Resultat 
dy O41 = DY Wi; (6) 


gefunden. Die linke Seite hat dabei (héchstens) 2n — 3 Glieder; die Anzahl in der rechten Seite 
ist gleich der Anzahl der Verbindungselemente 4;A,, d. h. héchstens 1/2 n (n — 1). 


vj? 
3. Fir ein raumliches System erhalt man analoge Resultate. Sind #, p Bi; und #7; die Winkel 


zwischen dem Vektor A?A? und der x-, y- und z-Achse, dann wird die potentielle Energie des 
Systems 


Say easel [(xj — x) cos B;; + (yj; — y;) cos Bi; + (z; — x) cos #;)?, (7) 


und die ersten drei Zeilen der Frequenzdeterminante 3n-ter Ordnung sehen so aus: 


> 6,1 cos? B,, — m, w? > 6; cos 0;, cos #5, >» G1 cos 0; cos 0; — Cy1 C08? Py 5 
Pa Ot 
2D) 1 cos J; cos 0,4 > 1 cos? 8, — m, wo? > ¢,, cos BH, cos P/, — cg, cos Pig COS Dig... 
© > c¢;, cos 0, cos. 0, > G1 cos 01 cos H, > ¢1 cos? B77, — m, w? — cg, cos Hy cos Pig... 


Teilt man diese Zeilen durch m,, die drei folgenden durch m, usw., betrachtet man wieder die 
Glieder der Hauptdiagonale, und beachtet, daB cos? J; ; + cos? 0;; + cos? Oj; = 1 ist, dann wird 
auch jetzt die Formel (4) gefunden. 


4. Die Elemente der Frequenzdeterminante hangen in einfacher Weise mit den Komponenten 
des Maxwellschen EinfluBtensors zweier Massenpunkte des Systems zusammen. Wenn man in der 
Gleichgewichtslage des (ebenen) Systems den Punkt A, eine Verschiebung mit den Komponenten 
x, ¥, gibt und die anderen Punkte samtlich an ihrer Stelle bleiben, dann wird in jedem Punkt 4, 
‘eine Kraft F,, hervorgerufen, deren Komponenten lineare Funktionen von x, und y, sind. Die 

‘Koeffizienten kann man aus (2) ablesen. Fiir s 4 | ist die Matrix ‘des Einflu8tensors 


3 : 
C), COS" @, C1, COS J}, SiN Y),. 


. s 2 
Ci, COS Y}, SIN Y), C;, SIN" Pis | 


(8) 


M,= | 
| 


Der Tensor ist nur dann ungleich null, wenn c,, 4 0, d. h. wenn A, und A, durch eine Feder ver- 

-bunden sind; er ist aber immer singular und entspricht der trivialen Tatsache, daf nur eine Ver- 
schiebung von A, in der Richtung A, A; eine Kraft in A, in derselben Richtung hervorruft. Fir 
s = | aber ist die Matrix 


— > 4; cos? 9; — Dd), COS Y; SIN Y}; 
i i 


Peek ee: 3 : (9 
Ii: Dd Gj C08 Yj; SID Yj — 2) Gj; sin? Hy; 


Die Schwarzsche Ungleichung sagt aus 


(Sa?) - (X bi) = (2 4,6), (10) 


20 


ase ee 


290 O. Bottema: Die Schwingungszahlen linear-elastischer Systeme Ingenieur-Archiy 


wobei das Gleichheitszeichen nur dann gilt, wenn die a; den b; proportional sind. Fir 
a, = cl? cosg,;, b= cf sing; 
folgt daraus, daB die Determinante von (9) nur dann gleich Null ist, wenn fiir ¢; 4 0 gilt, dab 
tg q,; fiir jedes i denselben Wert hat. Der Tensor (9) ist also nur dann singular, wenn samtliche 
mit A, verbundenen Punkte auf einer Geraden liegen. Ist das nicht der Fall, dann hat (9) zwei 
negative Eigenwerte, und es gibt dann in A, zwei senkrechte Hauptspannungsrichtungen. Ist a, 
der Winkel einer dieser Richtungen mit der x-Achse, dann ist zufolge einer bekannten Formel aus 
der analytischen Geometrie der Kegelschnitte 
Dd Ci Sin 2 jy 
te Qopes 2, 
tig, 2) Ci, COS 2 Pil qd) 
Nur wenn der Zahler und der Nenner in (11) beide gleich Null werden, sind die Hauptrichtungen 
unbestimmt und .4, ist ein ,,isotroper“‘ Punkt. 
Fiir den Raum gelten analoge Resultate. Der Tensor in A, ist 


| — > ¢; cos? #,, — > ¢; cos B,, cos i; — > ¢; cos B;, cos Be) 

=, ||| ua | 

M, = ne > 4; 008 Bj, cos 3,  — SG; cos? — >» aq; cos Bj, cos FB; : (12) 
— > a; cos Xf; cos 9;, — >) q; cos Hj, cos Bj1 — Sq, cos Of, | 


Er hat den Rang eins, wenn die mit A, verbundenen Punkte auf einer Geraden, und den Rang zwei, 
wenn sie in einer Ebene liegen. In jedem andern Fall hat er den Rang drei, und die Eigenwerte 
sind negativ. Sie bestimmen im allgemeinen in A, die drei Hauptrichtungen. 


Ay A; 


5. Wir wenden unsre Betrachtungen auf ein Beispiel an. Das System bestehe aus vier Punkten A, 
die in der Gleichgewichtslage ein Quadrat bilden. Jeder Punkt ist mit jedem anderen Punkt dace 
eine Feder verbunden. Die Punkte haben je die Masse m, die Feder je die Federkonstante c (Abb. 1) 

Wird die x-Richtung lings A,A4,, die y-Richtung langs A, A, angenommen, dann ist =e 


Ta=imy> (#+y5), 


2 
1 
fea 


1 
€ (x, — %9)? + = (x — %4)?-+ +e (ys = Ve)acae oe (¥1 — Ya)” 


1 1 
+ Ge (% 4% + ¥1 — Ys)? + 7 (ye — a — ya + 4)? 
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und fiir die Frequenzgleichung findet man, wenn man B=3— or setzt 
- C : 


—2 0 —l —1 0 0) 


1 
Pea oo ee ee 
0 
0 


7 sea) 2 ea 


0 Se Te) 6p a 1 
alee 0 £0 eye eh. —2 (isa ES) 
1 —l. 0-2 7I-S pw -.0 0| 
er U3 eee 0 3 
| (2 i See OT Ny 
oder 
(u — 3)? (u + I)*(u + 5) =0. -- (14) 


Wir haben also, wie zu erwarten war, eine dreifache Wurzel w? = 0 (welche die Bewegungen des 
starren Systems entspricht) und weiter als wesentliche Frequenzen eine vierfache Wurzel w? = w? 
= 2c/m und eine einfache Wurzel w? = w? = 4c/m. Die Amplituden der Schwingungen kénnen 
jetzt berechnet werden. Ist x,—= P,coswt, y; = Q; cos wt, dann ist fiir w? = w? der Rang der 
Frequenzdeterminante gleich vier, und es bestehen also vier linear unabhangige Gleichungen fiir 


P;und Q;. Wahlen wir P,, Q,, P, und Q, beliebig, dann hat man - 


Oe et Oy Pa O35 teem tee | (15) 
2Py= P, 0, Ee Qs » ‘ 2Q, = P,— Qj — P,—Q,. 
Fir @? = w} hat man 
P,=0,=Q,= Pyp=— P,=— P3=—Q=—Q; (16) 


und damit.ist die allgemeine Lésung des Schwingungsproblems gefunden. 
Wir bemerken noch, da8® immer 


PES eI OS EM : (17) 


ist, was der Tatsache entspricht, das der Schwerpunkt des Systems in Ruhe bleibt. 
Aus (11) geht hervor, da fiir jeden Massenpunkt die Hauptrichtungen mit den Diagonalen des 
Quadrats parallel sind. aes 


6. Die zweite Neubersche Formel sagt aus 


1 1 

Dr ee ake ee 
wo ,,00, die Eigenfrequenzen der Ersatzsysteme kennzeichnen, die sich durch Erstarrung aller Federn 
bis auf die jeweils fiir das Zustandekommen des Schwingungsvorganges mindest erforderlichen 
~ Federn bilden lassen‘. In unserem Beispiel kann man, damit ein schwingungsfahiges System tibrig 
bleibt, nicht mehr als vier Federn erstarren lassen. Es gibt also insgesamt 15 solche Ersatzsysteme, 
je von einem Freiheitsgrad. Wenn man aber einem Schwingungssystem Bedingungen auferlegt, 
dann wird immer die héchste Eigenfrequenz erniedrigt. Das heift: fiir die Frequenz @, des Ersatz- 
systems gilt @] << w? = 4c/m, und man hat also 


1 15 m 
a 


Das linke Glied von (18) hat aber den Wert 4 3 Se = a 


(Eingegangen am 18. August 1960.) 
Anschrift des Verfassers: Prof. Dr. O. Bottema, Delft (Holland), Julianalaan 132. 
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Berichtigung 


zu meiner Arbeit: 


Die Knickung der tordierten Welle mit Einzelkraft und kontinuierlicher Langs- 
kraft“ in Bd. 30, (1961), S.42 des Ingenieur-Archivs. 


Von H. Leipholz 


Seite 42, Uberschrift: kontinuierlicher (statt kontinuierlichen) ; 


Gleichung (4): t x ® = (statt t + &). 


Seite 43, Gleichung (5): in der dritten Gleichung = (state B) ; 


Seite 54, Zeile 2, letztes Glied in der letzten Klammer: 


2=r4s) Si, (stat — r,(s) S32) : 
Seite 56, Zeile 23: 


Fski 1 
0 : F 
f S37, 7, sin w s ds (sat f S2,7r, sin ws i) : 
0 0 


(Eingegangen am 29. Marz 1961). 
Anschrift des Verfassers: Dr. Horst Leipholz, Stuttgart-Untertiirkheim, Fiechtnerstr. 51. 
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Vor kurzem erschien: 


Praktische Mathematik 


as fiir Ingenieure und Physiker 


. Von Dr.-Ing. R. Zurmint, Darmstadt 


Dritte, verbesserte Auflage 


Mit 132 Abbildungen. XVI, 548 Seiten Gr.-8°. 1961. Ganzleinen DM 29,40 
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